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Introduction 


Introduction 


The three books (Parker 2011 and Parker 2012), Introductory, Intermediate and Advanced Maths for 
Chemists are tutorial workbooks intended for first year undergraduates taking a degree in Chemistry or a 
Chemistry based subject such as Chemical Engineering , Chemical Physics, Molecular Biology, Biochemistry 
or Biology. The texts may also be very useful for final year school or college students prior to them 
starting an undergraduate degree and their Chemistry teachers. The books are split into three in order 
to reduce file size and make handling on a laptop or tablet computer easier. Introductory Maths for 
Chemists roughly covers the first 8 weeks of semester 1; Intermediate Maths for Chemists the remainder 
of semester 1 and the beginning of semester 2; and Advanced Maths for Chemists the rest of semester 2. 
They each have chapter heading such as Week 1, Week 2 and so on. This is purely to help you self-pace 
your work on a weekly basis although I realize that the week numbers in Intermediate and Advanced 
Maths for Chemists will not correspond to your semester week numbers. 

This workbook is not a textbook! Use your Chemistry and Maths textbooks to find out the details about 
the areas covered. People will not really understand any subject, including Maths, until they can use it in 
a flexible way. Tutorials are a way of allowing you to practice your skills, in this case the Maths required 
by Chemists, so that the Maths will become easier with practice. 

As far as possible the workbook is organized on a weekly basis. Go through the examples and work 
out the solution yourself on paper then check the solution. A full solution is there for you to check that 
you are correct and to show my method of solving the problem. To begin with, the solutions give every 
single step but as you progress through the workbook the explanations become less detailed. When you 
do finally cover the Chemistry involved in the examples during your Chemistry degree you wont be 
blinded or scared by the Maths, as by then you will be happy playing around with equations. 

Advanced Maths for Chemists may be used with any Maths textbook, however, it is designed to interface 
with the textbook (Stroud and Booth 2007). Despite its name of Engineering Mathematics Stroud and 
Booths book covers all the Maths needed by all the sciences and engineering subjects. Advanced Maths 
for Chemists gives chemical examples of the Maths concepts. If you want to look up any first year 
chemistry then any General Chemistry textbook is useful but the textbook (Blackman, Bottle, Schmid, 
Mocerino & Wille 2012) is excellent. For later on in your course then (Atkins & de Paula 2009) and 
(Levine 2009) have many examples of the interplay between Chemistry and Maths mainly in the area 
of Physical Chemistry. 
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One final comment. A common mistake of many students is thinking you need to memorize all the 
equations you come across in any area of the subject. This is impossible but there are a very small number 
of equations that become familiar by usage and which you remember without really having to try. All the 
rest comes from being able to apply your Maths to these familiar equations (or to the equations supplied 
in an exam or from a textbook) and this enables you to get to your target. 
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1. Week 1: Chemical Data: Linear 
Least Squares Curve Fitting 

1.1 Least Squares Curve Fitting 

Least square fitting is usually made out to be more difficult than it really is or it is treated as a spreadsheet 
“black art”. I am going to spend some time explaining the fundamentals which are very easy to follow. 
Typically in a lab experiment we have 3 to 6 data points (x, y ) and we would like the best straight line 
through them. An excellent explanation of least squares fitting is given by Stroud and Booth 2007, p750, 
also you may want to look again at section 1.1.4 of the Introductory Maths for Chemists (Parker 2011) 
on linear graphs. As well as the linear function, 


y = mc+ c 


any function/(v.) that is capable of being plotted as a linear graph may have the data points fitted by 
the least squares method; some of the more chemically useful are listed (a and b are constants). A log 
function gives a straight line from a semi-log plot of y against ln(v) with gradient = a and intercept = b. 

y — a In (x)+Z? 

A power function gives a straight line from a log-log plot of In (y) against ln(x) with gradient = b and 
intercept = ln(a). 


_ b 

y - ax 

An exponential function gives a straight line from a semi-log plot of In (y) against v with gradient = 
In (b) and intercept = ln(a). 


y = ab x 

A reciprocal function gives a straight line from a semi-reciprocal plot of y against llx with gradient = 
b and intercept = a. 

_ b 

y - a + — 
x 

In Week 1 we cover linear curve fitting. Week 2 we will discuss non-linear curve fitting and fitting your 
own model functions to experimental data. 
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1.1.1 Linear Least Squares (LLSQ) Curve Fitting 

For a linear or straight line function we want to find the function, 

/ (x ) = m x + c 


which most closely fits our data points. That is, for a linear graph we need to discover the gradient m 
and the intercept c which gives this best fit to our data. Firstly we assume our values of the independent 
variable x have no error (or negligible error) but that the dependent variable y may have an appreciable 
random scatter. 


An individual data point (x, y.) has a residual, R , which is defined as the difference between the 
experimental datay. and the value ofy. that is predicted by the model equation straight liney. = mx. + c. 

R= y i ■ - {mx i + c) 


Removing the bracket the residual is, 


R - y . - m x l ~ c 


We take the square of the individual residuals, SR, to give equal weight to those data points which lie 
above the predicted line and those that lie below the line predicted of the model equation. 


SR= (y, - rax, - cj 


We take the sum the square of the residuals, SSR, of each of the data points to give an overall picture of 
the scatter of all the n data points. The sum of square of all the residuals SSR is, 

n 

SSR - ^ (y, - m x t - c J 

i= 1 

where n = number of data points. We will obtain the best fit for a straight line graph by minimizing the 
value of the sum of squares of the residuals by varying the values of x and c, hence the common name 
of a “least square” fit. This best fit line is called the linear least squares (LLSQ) line. The easiest way of 
curve fitting of equations to experimental data is to use a spreadsheet. You can use LibreOffice Calc, or 
OpenOfhce Calc, which are both open-source and free, or Microsoft Excel, see the References. In Calc 
or Excel you use an array function called LINEST (“lin est”) to fit a linear equation. 
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The LINEST algorithm works by an iterative process. It increases m by a small amount and then decreases 
m by the same small amount. It then compares the original SSR with these two altered values and chooses 
the value for m which leads to the smallest of these three SSR values. Holding this value of m constant 
for the moment, the algorithm then alters the intercept c by increasing and decreasing it and choosing 
the value of c which gives the smallest of the resulting three SSR values. LINEST then keeps repeating the 
above pair of cycles for the two variables until the final value of SSR is only changing by an insignificant 
amount. Let us look at a worked example as the best way of seeing how LLSQ curve fitting works. 


1.1.2 Spectrophotometer Absorption Measurement of Protein Concentration 



Figure 1.1: schematic diagram of a spectrophotometer. 
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A protein has been complexed with a molecule to form a compound that absorbs in the UV-visible 
range of the spectrum. In the Chemistry lab we use a spectrophotometer, see Fig. 1.1, to measure the 
absorption of the light. White light from a radiation source is selected for the wavelength at which the 
complex has a strong absorption. The selected wavelength has an intensity I Q which is passed through a 
known concentration of the complex, c mol L -1 , and a known length of the sample, l cm. The intensity 
of the light will be reduced by absorption of the complex to an intensity I. The log(/ 0 /I) is called the 
absorbance, A, of the solution and from the Beer-Lambert law, see below, the absorbance is proportional to 
the concentration of the absorbing molecule, c, the length of the sample, /, and the absorption coefficient 
of the molecule, e. The absorption coefficient of the molecule e varies with the wavelength of light and 
hence the need for using monochromatic light. 


A = log 



eel 


Fig. 1.2 is the initial spreadsheet of different known concentrations of the protein-complex and their 
absorbances. Fitting a line through the data points gives us a calibration curve for the complex. Obtaining 
the calibration curve means that we can measure an unknown concentration of the complex, e.g. obtained 
from a patient in a hospital. 



A 

B 

1 

Protein/pg 

A 

2 

0 

-0.0020 

3 

5 

0.0874 

4 

10 

0.1760 

5 

15 

0.2467 

6 

20 

0.3274 


Figure 1.2: initial spreadsheet for Beer-Lambert calibration line. 


To use the function LINEST we firstly select a blank area of the spreadsheet of 2 columns by 5 rows, lets 
assume this array of blank cells is D3 to E7. While this array of blank cells is still selected (high lighted) 
you enter your LINEST array-equation in the function window of the spreadsheet as, 


= LINEST (B2:B6 , A2:A6,1,1) 
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These four parameters within LINEST have the following meanings. The cells B2:B6 are the y-data and 
the cells A2:A6 are the x-data. The third parameter is a switch; if for good chemical reasons you want 
the line to go through the origin (0, 0) then set this parameter to 0, otherwise to find the intercept set 
it to 1. The fourth parameter is also a switch, if it is set to 0 then only the gradient and intercept are 
calculated, we will set it to 1 to obtain all the other statistics as well. These four parameters are separated 
by commas. Do not click on the arrow next to the function window and do not press return as this 
LINEST formula is an array formula, the array is from A2 to E7, instead enter the array formula with 
shift-control-return on a PC or command-return on a Mac. The blank cells you selected will now fill up 
with the LLSQ fitted parameters as in Fig. 1.3. 



A 

B 

C 

D E 

F 

1 

Protein/pg 

A 


LINEST Output 


2 

0 

-0.0020 


Gradient 

Intercept 


3 

5 

0.0874 

Parameter 

0.016362 

0.003480 


4 

10 

0.1760 

Std Dev 

0.000416 

0.005091 


5 

15 

0.2467 

r 2 

0.998067 

0.006573 

Std Dev of y 

6 

20 

0.3274 

F 

1549.289 

3 

Deg Freedom 

7 



SSreg 

0.066929 

0.000130 

SSR 


Figure 1.3: final spreadsheet of Beer-Lambert calibration. 


I have typed labels into the spreadsheet in columns C and F to help us understand the statistics produced. 
In cell D3 is the gradient m and E3 the intercept c. D4 is the standard deviation of the gradient o m and E4 
the standard deviation of the intercept a. D5 is the correlation coefficient r 2 which for a good fit should 
be a number close to but less than unity. Of the other statistics, the ones of interest are E6 the number of 
degrees of freedom which is the number of data points, 5 in this example, minus the number of coefficients 
to be fitted, 2 for m and c, so in this example the number of degrees of freedom is 3. E7 is the minimized 
value of the SSR of the fitted line which has been minimized by the LINEST algorithm of the spreadsheet. 

We quote the equation for our least squares line with the first decimal figure of the standard deviations 
(rounded if necessary) giving us the last decimal figure of the gradient and the intercept (also rounded 
if necessary). The LINEST fitted Beer-Lambert law calibration line is, 


A= (0.0164± 0.0004 |J g‘')(c[Jg)+ (0.003 ± 0.005) 


From our model equation (the Beer-Lambert law) the absorbance is dimensionless A = log(I 0 /I) and 
should pass through the origin. The graph of A against c is shown in Fig. 1.4 and the gradient m should 
have a gradient of units 1/c or pg _1 . The intercept has the units of A and is dimensionless. A graph should 
always be plotted to check that no mistakes have been made and there is not a rogue data point that 
needs to be either measured again or not plotted as it distorts the data too much. The line shown in Fig. 
1.4 is our LLSQ line from the LINEST spreadsheet. 
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Figure 1.4: Beer-Lambert law LLSQ fit. 


The relative precision of the gradient and the intercept are, 


m 


0.0004 

0.0164 


= 0.02439 = 2.4% 


°c 


c 


0.005 

0.003 


= 1.666 = 167% 


Clearly interpolating within the range 0 to 20 pg of protein for the unknown sample is valid, but 
extrapolating for an amount of protein greater than 20 pg should definitely not be attempted due to the 
large uncertainty in the intercept. Looking at the graph in Fig. 1.4 we might be tempted to try extrapolation 
but the size of the standard deviation of the intercept tells us that we would need to make experimental 
measurements at higher concentrations and not extrapolate the graph. 
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You may also wonder about forcing the LINEST line through the origin of the graph. This accomplished by 
selecting the cells where the statistics are displayed, D3:E7, altering the third parameter from 1 to 0 as below, 


= LINEST (B2:B6 , A2:A6,0, l) 


and then pressing control-shift-return on a PC. The spreadsheet now looks like Fig. 1.5. 



A 

B 

C 

D E 

F 

1 

Protein/pg 

A 


LINEST Output 


2 

0 

-0.0020 


Gradient 

Intercept 


3 

5 

0.0874 

Parameter 

0.016594 

0 


4 

10 

0.1760 

Std Dev 

0.000223 

#N/A 


5 

15 

0.2467 

r 2 

0.999275 

0.006119 

Std Dev of y 

6 

20 

0.3274 

F 

5515.195 

4 

Deg Freedom 

7 



SSreg 

0.206521 

0.000150 

SSR 


Figure 1.5: LLSQ fit for the line when forced through the origin. 


As we have now forced the line though the origin there is zero error in the intercept, and the standard 
deviation of the gradient is about half its previous size at 1.2%. 


a 


m 


m 


0.0002 

0.0166 


= 0.01205 = 1.2% 


The degrees of freedom has increased from 3 to 4 as the intercept is now a constant and not a variable. 
Perhaps the only thing to give us any slight concern is that now SSR ~ 0.00015 whereas previously it 
was ~ 0.00013. Despite the low standard deviation, extrapolation is still dangerous without any evidence 
that the curve continues to be linear, in fact further experiments show it starts to become slightly curved 
above about 20 pg of protein. The Beer-Lambert Law fit for our data when we force the line through 
the origin is, 


A= (0.01661 0.0002 |jg _1 )(c|Jg) 


Let us use LINEST to carry out some LLSQ curve fitting questions. 

1.2 Week 1 Tutorial Questions 

1.2.1 Question 1: Enthalpy of Fusion of n-Carboxylic Acids 

Fig. 1.6 shows the enthalpy of fusion AiT fusion of four linear saturated carboxylic acids as a function of n 
the number of carbon atoms in the molecule. 
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A 

B 

1 

n 

A//fusion / kJ mol 1 

2 

6 

15.06 

3 

8 

21.38 

4 

10 

29.40 

5 

12 

36.70 


Figure 1.6: enthalpy of fusion of n-carboxylic acids. 


Use LINEST to find the equation that correlates the enthalpy of fusion of the carboxylic acids with the 
number of C-atoms in their length and also find the precisions of the gradient and intercept. Organic 
chemists use this type of equation to correlate large tables of experimental data (not just the 4 pieces of 
data shown here). Finally draw a graph of the data with your LLSQ line shown. 

Tump to Solution 1 (see page 23) 


1.2.2 Question 2: Electrochemistry and Thermodynamics 
The following redox reaction, 

Zn + 2Fe(CN)j*~ (aq) -^Zn 2+ (aq) + 2Fe(CN)*~ (aq) 


was run as an electrochemical cell, Fig. 1.7. 

voltmeter 

Zn electrode 
Zn 2+ ions — 


.NaCI-agar gel 
salt bridge 


l/Pt electrode 
Fe(CN) 6 4_ ions 


Figure 1.7: schematic diagram of the electrochemical cell. 


A student measured the cell voltage E to the nearest millivolt at various temperature T, Fig. 1.8. 
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A 

B 

1 

77 K 

E/V 

2 

275 

1.255 

3 

277 

1.249 

4 

279 

1.244 

5 

281 

1.239 

6 

283 

1.234 

7 

285 

1.231 

8 

287 

1.228 

9 

289 

1.223 

10 

291 

1.219 

11 

293 

1.215 

12 

295 

1.211 

13 

297 

1.206 

14 

299 

1.202 

15 

301 

1.198 

16 

303 

1.193 

17 

305 

1.189 

18 

307 

1.185 

19 

309 

1.181 

20 

311 

1.177 

21 

313 

1.172 

22 

315 

1.169 

23 

317 

1.165 

24 

319 

1.16 

25 

321 

1.154 

26 

323 

1.15 

27 

325 

1.147 


Figure 1.8: cell voltage and temperature. 


The cell voltage is related to the entropy and enthalpy of the redox reaction by the equation, 


^ AH TAS 
nF nF 

Where n is the number of electrons transferred in the redox reaction [n- 2 in this case) and the Faraday 
constant F = 9.6485xl0 4 C mol' 1 . Type up your own spreadsheet and use LINEST to determine AH and AS 
from a linear graph of E against T. What are the precisions of AH and AS and also their relative precisions? 

Tump to Solution 2 (see page 24) 
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1.2.3 Question 3: The Visible Spectrum of Iodine 

The I 2 molecule absorb at the red-end of the visible spectrum and thus appears bluish-purple. Fig. 1.9 
shows the visible absorption spectrum taken using the teaching lab UV-visible spectrophotometer in 
Heriot-Watt University Chemistry Department. 



wavelength / nm 


Figure 1.9: visible spectrum of iodine. 


The fine structure is called a vibrational progression and it runs from at least 650 nm down to about 
500 nm when it becomes continuous i.e. non-quantized. The vibrational level v' around 500 nm at 
which it becomes non-quantized is at the dissociation limit of the excited electronic state and thus the 
absorption spectrum changes from a quantized spectrum to a continuous spectrum with increasing 
energy of the photons absorbed. This allows us to measure in the lab the energy to jump from v = 0 of 
the ground electronic state to the v' max the dissociation limit as shown in Fig. 1.10. Spectroscopy will be 
covered fully in your Chemistry lectures later in your degree. 



excited electronic state 
Franck-Condon region 


ground electronic state 
bond distance 


Figure 1.10: schematic PE curves for a diatomic molecule. 


In Fig. 1.11 the spreadsheet of the absorption spectrum has column B with the wavelengths A/nm of the 
peaks for the transitions from v = 0 to the various v' (column A) in the excited state. Wavelength is not 
proportional to energy but wavenumber G(v')/cm _1 is proportional to energy and so the peaks of the 
vibrational progression are given in wavenumbers in column C of Fig. 1.11. 
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A 

B 

C 

D 

E 

1 

v' 

/Jnm 

G(v')/cm 1 

AG (v'+'A) 

v'+y 2 

2 

16 

577.1 

17,328 



3 

17 

573.7 

17,431 

103 

16.5 

4 

18 

570.7 

17,522 

92 

17.5 

5 

19 

567.4 

17,624 

102 

18.5 

6 

20 

564.5 

17,715 

91 

19.5 

7 

21 

561.5 

17,809 

95 

20.5 

8 

22 

558.7 

17,899 

89 

21.5 

9 

23 

556.1 

17,982 

84 

22.5 

10 

24 

553.3 

18,073 

91 

23.5 

11 

25 

551.0 

18,149 

75 

24.5 

12 

26 

548.3 

18,238 

89 

25.5 

13 

27 

545.8 

18,322 

84 

26.5 

14 

28 

543.6 

18,396 

74 

27.5 

15 

29 

541.4 

18,471 

75 

28.5 

16 

30 

539.1 

18,549 

79 

29.5 

17 

31 

537.0 

18,622 

73 

30.5 

18 

32 

535.0 

18,692 

70 

31.5 

19 

33 

532.9 

18,765 

74 

32.5 

20 

34 

531.1 

18,829 

64 

33.5 

21 

35 

529.4 

18,889 

60 

34.5 

22 

36 

527.5 

18,957 

68 

35.5 

23 

37 

525.8 

19,019 

61 

36.5 

24 

38 

524.2 

19,077 

58 

37.5 

25 

39 

522.8 

19,128 

51 

38.5 

26 

40 

521.4 

19,179 

51 

39.5 

27 

43 

517.2 

19,335 



28 

44 

516.0 

19,380 

45 

43.5 

29 

45 

514.8 

19,425 

45 

44.5 

30 

46 

513.6 

19,470 

45 

45.5 

31 

47 

512.5 

19,512 

42 

46.5 

32 

48 

511.4 

19,554 

42 

47.5 

33 

49 

510.5 

19,589 

34 

48.5 

34 

50 

509.6 

19,623 

35 

49.5 

35 

51 

508.8 

19,654 

31 

50.5 

36 

52 

508.0 

19,685 

31 

51.5 

37 

53 

507.2 

19,716 

31 

52.5 

38 

54 

506.5 

19,743 

27 

53.5 

38 

55 

505.8 

19,771 

27 

54.5 

40 

56 

505.0 

19,802 

31 

55.5 

41 

57 

504.5 

19,822 

20 

56.5 

42 

58 

503.9 

19,845 

24 

57.5 

43 

59 

503.4 

19,865 

20 

58.5 

44 

60 

502.9 

19,885 

20 

59.5 

45 

61 

502.4 

19,904 

20 

60.5 

46 

62 

502.0 

19,920 

16 

61.5 


Figure 1.11: initial spreadsheet of iodine visible absorption spectrum. 
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Advanced Maths for Chemists 


Week 1: Chemical Data: Linear Least Squares Curve Fitting 


Type up your own spreadsheet for Columns A to E. The spectrum in Fig. 1.9 shows that we only have 
peaks for the middle part of the progression, we don’t know the wavelength of the v = 0 to V = 0 transition 
as it is beyond the range of our spectrometer and so we need to look at the energy differences between 
neighbouring peaks. Column D of your spreadsheet should calculate A G(v’+Vz) 

AG(v'+ 1 / 2 )= G(v'+l)- G(v') 

You will notice that there is missing data for transitions to V = 41 and 42. Such missing data is common in 
science, so think carefully about how your spreadsheet can best handle the missing data points for the next 
step. The summation of these AG(v'+Vz) is equal to the dissociation energy in the excited electronic state, Df 


D 0 ' = ^AG(v'+y 2 ) 


The best way of proceeding when you don’t have the data from outside the spectrometer’s range is to 
fit a straight line through the graph of A G(v+Vz) against (v'+Vz). From the equation for the line solve 
it to find the where the line crosses the v-axes i.e. the x-intercept. This extrapolation of a vibrational 
progression in spectroscopy is called the Birge-Sponer extrapolation. Extrapolation is always a worrying 
thing so finding the standard deviations of the gradient and intercept is very important. Once we have 
the x-intercept from the equation as well as the gradient and the y-intercept we can evaluate the D Q ' 
graphically as the area under the triangle of the AG(v’+V4) against (v'+Vz) line. 

Tump to Solution 3 (see page 26) 


1.2.4 Question 4: Arrhenius Plot 

Fig. 1.12 shows the experimental data from a student lab experiment to measure the first-order rate 
constant of a chemical reaction against temperature. 



A 

B 

1 

77K 

k / s 1 

2 

293 

0.00526 

3 

298 

0.00756 

4 

303 

0.0144 

5 

308 

0.0227 

6 

313 

0.0483 

7 

318 

0.0621 

8 

323 

0.119 


Figure 1.12: first-order rate constant and temperature data. 
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Advanced Maths for Chemists 


Week 1: Chemical Data: Linear Least Squares Curve Fitting 


The Arrhenius equation can be transformed into a linear equation by taking logs on both sides and 
plotting the linear equation In (A) against 1/T, this transformed graph is called an Arrhenius plot. 



E a \ 

/ \ / \ E „ 

k — Aex p 


InU) = InU) 


The intercept of the Arrhenius plot is c- ln(A) and the gradient m = -EJR. Type up your own spreadsheet, 
transform the data and find the pre-exponential factor A and the activation energy £ a and their precisions. 
Do you notice anything worrying about your LLSQ fitted line when you graph it? 

lump to Solution 4 (see page 29) 
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Advanced Maths for Chemists 


Week 1: Chemical Data: Linear Least Squares Curve Fitting 


1.3 Week 1 Tutorial Solutions 
1.3.1 Solution 1: Enthalpy of Fusion of n-Carboxylic Acids 
Fig. 1.13 shows the final spreadsheet using LINEST. 



A 

B 

C 

D 

E 

1 

n 

A//fusion / kJ mol -1 




2 

6 

15.06 


3.647 

-7.188 

3 

8 

21.38 


0.115 

1.070 

4 

10 

29.40 


0.998 

0.516 

5 

12 

36.70 


998.319 

2 

6 




266.012 

0.533 


Figure 1.1 3: enthalpy of n-carboxylic acids fitted by LINEST. 


Taking account of the standard deviations, the fitted equation for the enthalpy variation is, 


A// ihsion = (3.6±0.1 kJ mol” 1 )n - (7± 1 kJmol" 1 ) 


Notice how the intercept is much less precise (14% relative precision) compared to the gradient (3%) so 
extrapolation is not recommended. However, with more data (such as from a large tabulation of data) 
the precision of the LLSQ correlation line would improve markedly and become really useful in organic 
chemistry Also with more data you would find that you need two slightly different linear equations 
for n with even and odd numbers of C-atoms. Finally we should always check our results by plotting a 
graph to see if any of the points are rogue ones and to confirm that the fitted line is sensible, Fig. 1.14. 



Return to Question 1 (see page 16) 
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Advanced Maths for Chemists 


Week 1: Chemical Data: Linear Least Squares Curve Fitting 


1.3.2 Solution 2: Electrochemistry and Thermodynamics 



A 

B 

c 

D 

1 

77 K 

E/V 



2 

275 

1.255 



3 

277 

1.249 

-0.002121 

1.836256 

4 

279 

1.244 

0.000011 

0.003326 

5 

281 

1.239 

0.999347 

0.000847 

6 

283 

1.234 

36705 

24 

7 

285 

1.231 

0.026326 

0.000017 

8 

287 

1.228 



9 

289 

1.223 



10 

291 

1.219 



11 

293 

1.215 



12 

295 

1.211 



13 

297 

1.206 



14 

299 

1.202 



15 

301 

1.198 



16 

303 

1.193 



17 

305 

1.189 



18 

307 

1.185 



19 

309 

1.181 



20 

311 

1.177 



21 

313 

1.172 



22 

315 

1.169 



23 

317 

1.165 



24 

319 

1.16 



25 

321 

1.154 



26 

323 

1.15 



27 

325 

1.147 




Figure 1.15: LLSQ LINEST fit of cell voltage and temperature. 


The LLSQ line through the E versus T data with units is, 


E= (-0.002121 0.00001 VK^Xr K)+(1.8361 0.003 V) 


Plotting the graph lets see if we have made any miscalculation and whether there are any rogue points 
which may well need remeasuring in the lab, also a graph is required to complete your lab report, Fig. 1.16. 



Figure 1.16: LLSQ line for E versus T data. 
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Advanced Maths for Chemists 


Week 1: Chemical Data: Linear Least Squares Curve Fitting 


From the intercept for the graph and our thermodynamic equation we have, 

^ AH TAS 
nF nF 

AH = (- 1.836 ± 0.003 V)(2x 9.6485* 10 4 C mol' 1 ) 

AH = (-354.293, 4.-354.872, t-353.714) kJmol' 1 

We have used the fact that 1 coulomb-volt equals 1 joule. I have used an up and a down arrow to show 
the plus and minus values for A H at this intermediate stage and to reduce rounding errors an excess of 
significant figures is quoted. It is in the first decimal place, after rounding, that the three values vary thus, 

AH = (-354.3 ± 0.6) kJ mol' 1 

The gradient leads to, 

AS = (-0.00212 ± 0.00001 V K' ')(2* 9.6485* 10 4 Cmol' 1 ) 

AS= (-409.0964, |-411.0261, | “ 407.1667)JK" 1 mol' 1 
AS = (- 409 ± 2) J K 1 mor 1 

It is noticeable that the relative standard deviation of AH is -0.2%, whereas AS is -0.5%. We have tried 
to reduce rounding errors in the different stages of the data processing and to show the importance of 
finding the standard deviations of the gradient and intercept and not just m and c themselves. 

Return to Question 2 (see page 17) 
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Week 1: Chemical Data: Linear Least Squares Curve Fitting 


1.3.3 Solution 3: The UV-Visible Spectrum of Iodine 



A 

B 

C 

D 

E 

F 

G 

H 

I 

1 

v' 

2/nm 

G(v' )/cm 1 

AOKv'+Vi) 

v'+y 2 

AGfv’+Vi) 

V'+‘/2 



2 

16 

577.1 

17,328 







3 

17 

573.7 

17,431 

103 

16.5 

103 

16.5 

-1.893 

130.174 

4 

18 

570.7 

17,522 

92 

17.5 

92 

17.5 

0.043 

1.750 

5 

19 

567.4 

17,624 

102 

18.5 

102 

18.5 

0.980 

3.822 

6 

20 

564.5 

17,715 

91 

19.5 

91 

19.5 

1983.872 

41 

7 

21 

561.5 

17,809 

95 

20.5 

95 

20.5 

28983.986 

599.002 

8 

22 

558.7 

17,899 

89 

21.5 

89 

21.5 



9 

23 

556.1 

17,982 

84 

22.5 

84 

22.5 



10 

24 

553.3 

18,073 

91 

23.5 

91 

23.5 



11 

25 

551.0 

18,149 

75 

24.5 

75 

24.5 



12 

26 

548.3 

18,238 

89 

25.5 

89 

25.5 



13 

27 

545.8 

18,322 

84 

26.5 

84 

26.5 



14 

28 

543.6 

18,396 

74 

27.5 

74 

27.5 



15 

29 

541.4 

18,471 

75 

28.5 

75 

28.5 



16 

30 

539.1 

18,549 

79 

29.5 

79 

29.5 



17 

31 

537.0 

18,622 

73 

30.5 

73 

30.5 



18 

32 

535.0 

18,692 

70 

31.5 

70 

31.5 



19 

33 

532.9 

18,765 

74 

32.5 

74 

32.5 



20 

34 

531.1 

18,829 

64 

33.5 

64 

33.5 



21 

35 

529.4 

18,889 

60 

34.5 

60 

34.5 



22 

36 

527.5 

18,957 

68 

35.5 

68 

35.5 



23 

37 

525.8 

19,019 

61 

36.5 

61 

36.5 



24 

38 

524.2 

19,077 

58 

37.5 

58 

37.5 



25 

39 

522.8 

19,128 

51 

38.5 

51 

38.5 



26 

40 

521.4 

19,179 

51 

39.5 

51 

39.5 



27 

43 

517.2 

19,335 



45 

43.5 



28 

44 

516.0 

19,380 

45 

43.5 

45 

44.5 



29 

45 

514.8 

19,425 

45 

44.5 

45 

45.5 



30 

46 

513.6 

19,470 

45 

45.5 

42 

46.5 



31 

47 

512.5 

19,512 

42 

46.5 

42 

47.5 



32 

48 

511.4 

19,554 

42 

47.5 

34 

48.5 



33 

49 

510.5 

19,589 

34 

48.5 

35 

49.5 



34 

50 

509.6 

19,623 

35 

49.5 

31 

50.5 



35 

51 

508.8 

19,654 

31 

50.5 

31 

51.5 



36 

52 

508.0 

19,685 

31 

51.5 

31 

52.5 



37 

53 

507.2 

19,716 

31 

52.5 

27 

53.5 



38 

54 

506.5 

19,743 

27 

53.5 

27 

54.5 



38 

55 

505.8 

19,771 

27 

54.5 

31 

55.5 



40 

56 

505.0 

19,802 

31 

55.5 

20 

56.5 



41 

57 

504.5 

19,822 

20 

56.5 

24 

57.5 



42 

58 

503.9 

19,845 

24 

57.5 

20 

58.5 



43 

59 

503.4 

19,865 

20 

58.5 

20 

59.5 



44 

60 

502.9 

19,885 

20 

59.5 

20 

60.5 



45 

61 

502.4 

19,904 

20 

60.5 

16 

61.5 



46 

62 

502.0 

19,920 

16 

61.5 






Figure 1.1 7: LINEST fit for l 2 by Birge-Sponer extrapolation. 
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Week 1: Chemical Data: Linear Least Squares Curve Fitting 


Because of the missing data for V 41 and 42 the energy difference of 43-40 is anomalously large at 
155.7 cm -1 at D27 and has been deleted from the spreadsheet, otherwise it would distort the curve 
fitting. If we try to use LINEST directly on columns D and E then it throws up an error. This is because 
of the blank cells due to the missing data are treated as text and you cant have text in the middle of 
a numerical array. This is easily fixed. Columns D and E are selected and copied. Selecting the blank 
column F click on Edit then Paste Special, and the window shown in Fig. 1.18 appears in Calc (or the 
equivalent window in Excel). Deselect “Paste all” and select “Text”, “Numbers” and “Skip empty cells” 
and leave the other options alone. Press OK and the amended copy of D and E is copied into F and G 
but without the blank cells. 
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Advanced Maths for Chemists 


Week 1: Chemical Data: Linear Least Squares Curve Fitting 



Figure 1.18: paste special window in Calc 

Select a 2 column by 5 row blank area and enter the LINEST equation into the formula window while 
the cells are still selected and then press Shift-Control-Return on a PC or Command-Return on a Mac, 

= LINEST(F3:F45 , G3:G45,1 , l) 

Taking account of the standard deviations the equation of the Birge-Sponer extrapolation is thus, 

A G(v'+ 1 / 2 )= (- 1.89± 0.04cm’ 1 )(v'+V 2 ) + (l30± 2cm’ 1 ) 


The relative precisions are about 2% for the gradient and about 2% for the intercept, which is good 
precision for an undergraduate lab experiment and reflects the both scatter seen in the graph of Fig. 1.19 
but also reflects the advantage of the large number of data points in the experiment. 



Figure 1.19: LINEST fit for Birge-Sponer extrapolation for l 2 . 

The I 2 molecule dissociation limit of the excited electronic state, i.e. the intercept on the v-axis (rather 
than the y- intercept) is given by the equation of the line and is when AG(v' + Vi) = 0 we have, 

0 = (- 1.89± 0.04cm _1 )(v ' lirrit +1 /^) + (l30± 2cm' 1 ) 
v' linit +!/ 2 = 68.8 T 71.4 | 66.3cm -1 
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Where I use and up and down arrow to show the maximum and minimum values obtained. 

v 'i„m +1/2 “ 69± 3 cm” 1 

Where the precision of the v-intercept is found by firstly using the maximum value of m and the minimum 
value of c; and secondly the minimum value of m and the maximum value of c. We have treated (V+V 2 ) 
as a continuous variable as we must in order to use the least squares fitting procedure and to find the 
area under the curve. However as v' is actually quantized then the last quantized transition would be to 
v'limit = 69 ± 3 and any non-quantized transition above 69 leads to dissociation into two iodine atoms. 
Note that the scatter in our experiment is considerably limiting the precision but this might not be 
obvious without the use of LINEST. The dissociation energy of the excited state DJ is the area under 
the curve which is triangular and so, 

D 0 '= V 2 (69± 3 )(l30± 2cm _1 )= 4485 f 4752 f 4224cm" 1 
D 0 '= 4500± 300cm' 1 

As 1 cm -1 = 1.9864xl0“ 23 J and N A = 6.02214xl0 23 mol -1 the students lab experiment gives, 

zy = (4500± 300cm" 1 )(l.9864x 10" 23 jX6.02214x lO^mol" 1 ) 

D 0 '= (53.8 t 57.4 f 50.2)kJ mol" 1 
D 0 ' = (54± 4)kJ mol' 1 

The relative precision of our experimental excited state dissociation energy if about 7%. We have learnt 
how to handle missing experimental data in the middle and beginning of a data set. 

Return to Question 3 (see page 19) 


1.3.4 Solution 4: Arrhenius Plot 

Type up your own spreadsheet of the T, k, 1 IT and ln(k) in columns A, B, C and D. Use LINEST to give 
a spreadsheet similar to Fig. 1.20. 
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A 

B 

C 

D 

E 

F 

1 

77K 

k / s 1 

l/T 

ln(yt) 



2 

293 

0.00526 

0.003413 

-5.248 

-9986.79 

28.73 

3 

298 

0.00756 

0.003356 

-4.885 

411.32 

1.34 

4 

303 

0.0144 

0.003300 

-4.241 

0.99 

0.11 

5 

308 

0.0227 

0.003247 

-3.785 

589.52 

5 

6 

313 

0.0483 

0.003195 

-3.030 

7.79 

0.06609 

7 

318 

0.0621 

0.003145 

-2.779 



8 

323 

0.119 

0.003096 

-2.129 




Figure 1.20: LINEST fit for an Arrhenius plot. 


The LINEST fit to the plot of the log form of the Arrhenius equation is, 


In (A) = In (A)- 


RT 


In it = 29± 1 + — 10000±400K 


-K 


From the intercept the pre-exponential factor A is found by taking the anti-log, i.e. exp(29 ± 1). I have 
used some extra decimal places to reduce any rounding errors. As it is an exponential scale the error 
limits are not symmetrical around the mean value. The pre-exponential factor from our LLSQ fitted 


experimental data is, 
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A = (3.93 t 10.69 i 1.45) x lO'Y 1 


A= (4 +7 - 3) x 10 12 s _1 


The activation energy £ a is obtained from the gradient as, 

£ 

-~z L = -(10000± 400 K) 

R 

E a = (10000± 400K)(8.314J K~ ^oT ! ) 
E a = (83.14 t 86.47 f 79.8l)kJ mol" 1 


E a = (SStSjkJmol" 1 


The relative errors of the pre-exponential factor are between +175% and -75% and so extrapolation 
should be avoided. However, the relative error of the activation energy is fairly small at 4%. Fig. 1.21 
shows the Arrhenius plot and the LLSQ fitted equation in red. 



0.0031 0.0032 0.0033 0.0034 
Figure 1.21 : Arrhenius plot with LINEST fitted lines. 


The red line in Fig. 1.21 is the LINEST statistically more reliable line taking the standard deviations into 
account by plotting ln(k) = 29 - (10000 IT). The blue line is the LINEST fit without taking into account 
the standard deviations of the gradient and intercept by plotting ln(k) = 28.7 - (9987/T). The blue line 
is probably what you would obtain by using a plastic ruler and drawing a line by “eye” (which is how 
such plots were originally treated) and it is a pretty good fit! 

So what is going on? We are not plotting a y versus v graph but a transformed graph of ln(y) versus l/x. 
In Fig. 1.22, using the current data points, we have plotted the reciprocal of the temperature against the 
temperature which gives a linear plot with an excellent goodness of fit of r 2 = 0.9992 over this limited 
temperature range. 
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0.0034- 
0.0033 ■( 1 /7) /K -i 
0.0032- 
0.0031 


77K 


290 295 300 305 310 315 320 

Figure 1.22: 1 1T versus T for the Arrhenius experimental data. 


Conversely, Fig. 1.23 shows that the plot of the natural log of the rate constant against the rate constant 
and the graph is markedly curved for the current data. 



Figure 1.23: ln(/c) versus k for the Arrhenius experimental data. 

At large values of k then In (k) is nearly linear with /c, but a small values of k the log increases rapidly. These 
small rate constants will also be more difficult to measure accurately and will have a larger relative error. 
Perhaps transforming equations in order to obtain a linear Arrhenius plot may not be using the data to 
its best advantage. We need a non-linear fitting method, a NLLSQ method. This is what we look at next. 

Return to Question 4 (see page 21) 
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2 Week 2: Chemical Data: Non- 
Linear Least Squares Curve 
Fitting 

2.1 Nonlinear Least Squares (NLLSQ) Curve Fitting 

There are many equations and laws in Chemistry which may be plotted as linear graphs, as well as the 
obvious y = mx + c, e.g. log functions/(v) = aln(x) + b; power functions f(x) = ax b ; exponential functions, 
f(x) = a\f\ reciprocal functions,/(x) = a + b/x. The transformed functions can be treated by the LLSQ 
LINEST approach which was covered in Week 1. As we saw with the Arrhenius plot such transformations 
may not use the data to its best advantage. There are also many equations and laws that cannot be plotted 
as straight lines e.g. a double log function,/(v) = aln(v) + bln(x). 

In both LibreOffice Calc and OpenOffice Calc the algorithm Solver is part of the software but in Excel 
you may have to install Solver for NLLSQ, see the appropriate web-site. The simplest way of showing 
how to use this system is using a worked example. I will use LibreOffice Calc but anyone using Excel 
will find similar instructions. 
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2.1.1 NLLSQ Worked Example: Arrhenius Data 

We are going to use the data from Section 1.3.4 to directly plot k against T and fit a NLLSQ line using 
the Arrhenius equation as our model equation. Fig. 2.1 shows column A as our v variable and column B 
as y. Column C is the calculated value of y using the equation Arrhenius equation, k(c ale) = A exp(-£ / 
RT). So C2 contains the formula, 


= $G$3*EXP(- $G$4/($G$5*A2)) 



A 

B 

C 

D 

E 

F 

G 

1 

77 K 

k / s 1 

k{ calc) 

R 

SR 



2 

293 

0.00526 

0.999590 

-0.994330 

0.988691 



3 

298 

0.00756 

0.999596 

-0.992036 

0.984136 

A = 

1.00E+000 

4 

303 

0.0144 

0.999603 

-0.985203 

0.970625 

E = 

a 

1 

5 

308 

0.0227 

0.999610 

-0.976910 

0.954352 

R = 

8.314 

6 

313 

0.0483 

0.999616 

-0.951316 

0.905002 



7 

318 

0.0621 

0.999622 

-0.937522 

0.878947 



8 

323 

0.119 

0.999628 

-0.880628 

0.775505 



9 




SSR = 

6.457259 




Figure 2.1 : Arrhenius data with A = 1 and £ A = 1. 


In this function the arbitrary constants (i.e. the general name for the as yet unknown but soon to be 
determined constants) A and £ a , are entered as absolute cell references $G$3 and $G$4. We will return 
to the values these constants in a moment. The value for the gas constant R is also entered as cell $G$5. 

Column D are the individual residuals R , e.g. in cell D2 the following is entered, =B2-C2. The residuals 
squared SR are in column E e.g. cell E2, is =D2 A 2. Cells C2, D2 and E2 are copied and then pasted 
into the cells below. The sum of squares of the residuals SSR in E9 is =SUM(E2:E8). We will find the 
best straight line by minimizing the SSR by varying the unknown arbitrary constants A and E . This 
minimization of SSR occurs iteratively with A and £ a being altered slightly and randomly and seeing 
which change leads to a lower value for SSR. In the Solver algorithm we can change the total number of 
iterations allowed (called the learning cycles) and also the test of when we have reached the minimum 
for SSR (called the stagnation limit). 

There may be several minima in the three-dimensional space (the two arbitrary constants A and £ a and 
the SSR), but we want to find the lowest minimum, the so-called global minimum. Unlike the linear 
least squared method we need to specify the initial guesses for the two constants using our Chemistry 
knowledge so that we find the global minimum and not get stuck in a local minimum which may solve 
the Maths but is Chemical nonsense. 
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For the initial starting values of the arbitrary constants if the values of unity (1) are typed into G3 and 
G4 the SSR = 6.457 and so unity is a poor initial estimate for the constants. In this particular example we 
already have approximate values for the pre-exponential factor and the activation from the LLSQ LINEST 
fitting of the transformed equation. So lets use those values. Typing 1E12 into G3 and 80000 into G4 
gives us a much better value for SSR of 0.000155 as Fig. 2.2.. This is a better initial estimate than having 
them as unity. Now we can use the NLLSQ program Solver to automate this iterative process of altering 
A and £ a to find the global minimum pre-exponential factor and activation energy by a direct NLLSQ fit. 



A 

B 

C 

D 

E 

F 

G 

1 

77K 

k / s 1 

£(calc) 

R 

SR 



2 

293 

0.00526 

0.005463 

-0.000203 

4.1393E-08 



3 

298 

0.00756 

0.009479 

-0.001919 

3.6833E-06 

A = 

1.00E+012 

4 

303 

0.0144 

0.016150 

-0.001750 

3.0632E-06 

E = 

a 

80000 

5 

308 

0.0227 

0.027044 

-0.004344 

1.8870E-05 

R = 

8.314 

6 

313 

0.0483 

0.044546 

0.003754 

1.4091E-05 



7 

318 

0.0621 

0.072233 

-0.010133 

1.0268E-04 



8 

323 

0.119 

0.115388 

0.003612 

1.3047E-05 



9 




SSR = 

0.000155 




Figure 2.2: Arrhenius data with reasonable values for A and f A . 


Notice that the value of the gas constant is in J K _1 mol -1 and for the units to cancel in the exponential 
the activation energy must also be in the base unit of J mol -1 and is typed in as 80000. I am using 
LibreOffice Calc but Excel or OpenOffice have similar commands, place the cursor on the SSR cell i.e. 
E9. Click on Tools and then click on Solver. The Solver window that appears as in Fig. 2.3 with the 
Target cell shown as $E$9. 



Figure 2.3: Solver Window. 
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As we need to minimize the SSR we click on “Minimum” in the Solver window. Solver needs to be told 
where the cells are that it may vary in order to achieve the NLLSQ result. In the Solver window “By 
changing cells” use the drop down windows to select $G$3:$G$4. We can speed things up by telling 
Solver that the pre-exponential factor and activation energy are both positive values. As we know from 
our Chemistry that the Arrhenius parameters are both positive, using the “Limiting conditions” drop 
down windows set $G$3 >= 0 and $G$4 >= 0 and this will speed the iterative process. Finally we choose 
the algorithm to use, click on “Options” and select “DEPS”, then click on “Learning cycles” and edit it 
from 2,000 to 2,000,000 and then click on “Stagnation limit” and edit it from 70 to 70,000. Click OK 
and then click “Solve” in the Solver window. 


On my laptop it took 11.7 seconds to reach stagnation and used 1,374 cycles of iteration. The final 
spreadsheet is shown in Fig. 2.4. 



A 

B 

C 

D 

E 

F 

G 

1 

77 K 

k / s" 1 

k{ calc) 

R 

SR 



2 

293 

0.00526 

0.004734 

0.000526 

0.000000 



3 

298 

0.00756 

0.008427 

-0.000867 

0.000001 

A = 

4.00E+012 

4 

303 

0.0144 

0.014718 

-0.000318 

0.000000 

E = 

a 

83726 

5 

308 

0.0227 

0.025244 

-0.002544 

0.000006 

R = 

8.314 

6 

313 

0.0483 

0.042559 

0.005741 

0.000033 



7 

318 

0.0621 

0.070582 

-0.008482 

0.000072 



8 

323 

0.119 

0.115238 

0.003762 

0.000014 



9 




SSR = 

0.000127 




Figure 2.4: Arrhenius data with NLLSQ fitted Arrhenius parameters. 


SSR has dropped to 0.000127 and the NLLSQ values for the arbitrary constants are A = 4xl0 12 s 1 and 
£ a = 83.726 kj mol -1 . Remember to save the spreadsheet. 

Noting these values for our two arbitrary constants (A and EJ the question arises what are their precisions? 
You manually vary each one of the constants in turn and look at the subsequent variations in the SSR. 
This is easy with a spreadsheet. We increase the first digit in the first constant A and look at the change 
in SSR. If SSR alters by more than a factor of 10, press the undo icon and then move on to the next 
lower digit. When you find a digit that alters the SSR by less than a factor of ten you may take that as 
a reasonable measure of the precision of the NLLSQ fit for that constant (A). You repeat the procedure 
for the other constant E . Remember dont save this “working” spreadsheet. In this way you quickly find 
that A = (4 ± l)xl0 12 s -1 and E & = (83.7 ± 0.1) kj mol -1 . The Arrhenius equation obtained by the NLLSQ 
fitting of the rate constant and temperature data is shown below and the data is plotted in Fig. 2.5. 


k — (4±lXl0 12 s l ) exp 


—83700 ±100 J mol 1 
RT 
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Figure 2.5: rate constant and temperature with NLLSQ fitted curve. 

The raw data is best used with a NLLSQ approach rather than a transformation to a linear equation 
which distorts the data and the relative precisions of the data points. The precisions of the NLLSQ 
arbitrary constants are found in the same way as we have been using for the LLSQ precisions. The 
NLLSQ gives more precise (reliable) values for the Arrhenius parameters. LINEST (section 1.3.4) gave 
us A = (4 + 7 -3)xl0 12 s _1 and £ a = (83 ± 3) kj mol“ 1 for the same data. However, LINEST did give us 
useful initial starting points for the arbitrary constants in the iterative process. It is worth stressing that 
we don’t need to carry out LINEST before using Solver, but we can directly go to the NLLSQ method, 
as in the following questions. 

But in general we can now check any model equation against our data to see whether the equation does 
in fact fit the experimental data. 





In the past four years we have drilled 



81,000 km 



That's more than twice around the world. 
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to help customers find and produce oil and gas safely. 
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A ^ 
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2.2 Week 2 Tutorial Questions 
2.2.1 Question 1: Boiling Points of n-Alkanes 
A typical n-alkane (w-nonane) is shown in Fig. 2.6 



Figure 2.6: n- nonane. 


The boiling points (BP K) of some liquid n-alkanes of formula C n H 2n+2 are believed to fit the formula 
shown below. This is an example of the type of equation which organic Chemists use to model large 
amounts of tabulated data. Clearly here we are using a very small data sample but increasing the number 
of data points will increase the precision and usefulness of the equation and the model equation will 
then be a convenient summary. 


BP= Alog(n + B)~ C 


Type up your own spreadsheet and solve it for the arbitrary constants (A, B , and C) to find the model 
equation from the data in Fig. 2.7 and then interpolate to find the boiling point of n-nonane and the 
precision of the predicted boiling point. 



A 

B 

1 

n 

BP/ K 

2 

5 

309 

3 

6 

342 

4 

8 

399 

5 

10 

447 

6 

11 

469 

7 

12 

489 


Figure 2.7: boiling points for some n-alkanes. 


Tump to Solution 1 (see page 43) 
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2.2.2 Question 2: Heat Capacity of Solid Lead 

The thermodynamic properties of solids and liquids are accurately modelled by the Shomate equation. 
For solids the Shomate equation is, 


C p = A + Bt + Ct 2 + Dt 3 + pr 


Where C p is in J K _1 mol -1 and t=(T K)/1000. Fig. 2.8 shows the measured heat capacity C p of solid Pb 
in column C from temperatures near room temperature at 300 K up to 600 K near its boiling point in 
column A, and in column B the temperatures t. 



A 

B 

C 

1 

77 K 

t = 771000 

C(p)/J K 1 mol- 1 

2 

300 

0.30 

26.7 

3 

350 

0.35 

27.3 

4 

400 

0.40 

27.7 

5 

450 

0.45 

28.1 

6 

500 

0.50 

28.6 

7 

550 

0.55 

29.0 

8 

600 

0.60 

29.4 


Figure 2.8: measured C p for Pb at various temperatures. 


Set up your own spreadsheet to solve the NLLSQ fitting for the Shomate equation for solid lead. 


lump to Solution 2 (see page 45) 


2.2.3 Question 3: Ionization of He 


Electron 


Ion repeller 
beam 


Gas inlet 

I 



Heated 
tungsten wire 


-^Electron trap 
Ion beam 


Figure 2.9: mass spectrometer electron-ionization source. 
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Fig. 2.9 shows a stainless steel source of a mass spectrometer which is inside a vacuum chamber at a 
pressure of about 10“ 9 atm. A beam of electrons (from a heated tungsten filament) is passed through a 
low pressure of helium gas (~10“ 7 atm.) and the probability of ionization, called the cross section S (in 
angstrom squared A 2 = 10“ 20 m 2 ) is measured as a function of the colliding electrons energy, £, (measured 
in a unit called the electron volt, eV). 


The impacting electron losses kinetic energy in the collision from KE 1 to KE2 and the energy lost is 
required to both ionize the helium atom and give the liberated slow electron its kinetic energy of KE3. 
Helium ionizes at E = 24.59 eV and then S increases to a broad plateau at E = 60-80 eV after which 
S decreases slowly with E. There are good theoretical reasons for believing that our high energy data, 
beyond the plateau region, should fit the equation, 



Where A and B are unknown coefficients. From the data in Fig. 2.10 design a spreadsheet and use Solver 
to find the unknown coefficients and then plot the data with its NLLSQ line. Using your equation what 
is the cross section for ionization of helium with 485 eV electrons and what is the precision of this 
cross-section? 
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A 

B 

1 

E/eV 

s/A 2 

2 

150 

0.43 

3 

175 

0.40 

4 

200 

0.39 

5 

250 

0.37 

6 

300 

0.34 

7 

350 

0.31 

8 

400 

0.28 

9 

450 

0.27 

10 

500 

0.26 

11 

550 

0.24 

12 

600 

0.23 

13 

650 

0.21 

14 

700 

0.21 

15 

750 

0.19 

16 

800 

0.19 

17 

850 

0.17 

18 

900 

0.17 

19 

950 

0.17 

20 

1000 

0.16 


Figure 2.10: ionization cross sections for He. 


Tump to Solution 3 (see page 47) 


2.2.4 Question 4: Enzyme Reaction Kinetics 

The rates of enzyme reactions literally determine our health and well-being, it is extremely important 
both personally and industrially (the pharmaceutical industry). Adenosine triphosphate (ATP) is used by 
all cell to transfer energy which is accessed by hydrolysis with myosin. Fig. 2.11 shows the molecule ATP. 



Figure 2.11 : ATP N = blue, C = dark grey, H = grey, P = orange, 0 = red. 
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For a constant enzyme concentration, the rates of many enzyme reactions follow the Michaelis-Menten 
equation, 

V maxl S ] 

^M+[S] 

where v is the rate of reaction and [S] is the substrate molecules concentration, for this reaction the 
substrate is ATP. K u is the Michaelis constant, characteristic of the enzyme and substrate combination. 

At large concentrations of substrate when [S] » K u then the denominator (K M + [S]) ~ [S] and so the 
rate plateaus off to v ~ v max . When v = v m J2 then K u = [S] 1/2 with [S] % being the substrate concentration 
corresponding to v max /2. The Michaelis-Menten equation is a nonlinear equation but may be rearranged 
into a double reciprocal form called the Line weaver-Burk equation. 

I = 1 , 

v v max v M [S] 


Plotting 1/v against 1/[S] gives a straight line graph called the Lineweaver-Burk plot. This transformed 
linear equation is very useful in distinguishing different types of enzyme inhibition. However, transforming 
the equation into a linear form distorts the structure of the data and the relative errors as it is strongly 
biased in fitting the data at low values of [S] and low v which are the least precise part of the data set with 
the largest relative error, it is inferior to a NLLSQ fitting. Fig. 2.12 is the data obtained by a student using 
a constant myosin enzyme concentration for the hydrolysis of ATP at 25°C and pH = 7.0. The [ATP]/10“ 5 
mol L _1 and v/10 -6 mol L _1 s _1 however they are shown as [ATP] and v for clarity in the spreadsheet. 



A 

B 

1 

[ATP] 

V 

2 

7.5 

0.06 

3 

12.5 

0.10 

4 

25 

0.13 

5 

32.5 

0.15 

6 

50 

0.16 

7 

75 

0.18 

8 

100 

0.19 

9 

125 

0.19 

10 

150 

0.20 

11 

175 

0.21 

12 

200 

0.20 

13 

225 

0.21 

14 

250 

0.20 

15 

275 

0.21 

16 

300 

0.21 

17 

325 

0.21 

18 

350 

0.20 


Figure 2.12: [ATP] and its rate of hydrolysis by myosin. 
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Plot the raw data as a Michaelis-Menten equation (v against [S]) and use the fact that the curve goes to 
v max at large [S] and when v = v m J2 then K M = [S] . Use these values to help with your starting values 
of the two constants. Solve the data by a NLLSQ fitting of the Michaelis-Menten equation and find a 
measure of the precision of the constants. 

Tump to Solution 4 (see page 48) 


2.3 Week 2 Tutorial Solutions 

2.3.1 Solution 1: Boiling Points of n -Alkanes 



A 

B 

C 

D 

E 

F 

G 

1 

n 

BP/ K 

BP calc 

R 

SR 



2 

5 

309 

309.034 

-0.034 

0.001 

A = 

724.533 

3 

6 

342 

341.949 

0.051 

0.003 

B = 

4.069 

4 

8 

399 

398.960 

0.040 

0.002 

C = 

384.741 

5 

10 

447 

447.209 

-0.209 

0.044 



6 

11 

469 

468.815 

0.185 

0.034 



7 

12 

489 

489.033 

-0.033 

0.001 



8 








9 




SSR = 

0.0842 




Figure 2.13: n-alkane boiling points for NLLSQ fitting. 
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Fig. 2.13 shows the spreadsheet for using Solver, column C has the calculated boiling point using the 
model equation with C2 having the equation, 

= $G$2* LOG(A2 + $G$3) - $G$4 

typed in the formula window. Cell D2 has =B2 - C2, and E2 has = D2 A 2 entered. The rests of the columns 
C, D and E are filled in by copying and pasting. The SSR is in cell E9 with =SUM(E2:E7) entered. 

We need to find a reasonable initial starting values for the arbitrary constants. If you type in unity (1) 
in G2, G3, G4 then the SSR is over a million and is clearly wrong. The boiling points are in the range 
of 300-400 K which arises very approximately as the difference between A and C. So typing in a few 
values by hand then A = 700 and C = 300 gives SSR = 8234, which is a lot smaller than a million. The 
constant B fine tunes the boiling point for changes in chain length of the alkane. So stating from B = 1 
try some numbers and you find that B - 2 gives an SSR = 765. This should be a good enough starting 
point for the DEPS algorithm to find the global minimum. 

Placing the cursor on cell E9 click on Tools, Solver. The target cell should be $E$9, click on minimize, 
select cells G2 to G4 to change. In the Options of the Solver window make sure the algorithm used is 
DEPS, edit learning cycles from 2,000 to 2,000,000 and the stagnation limit from 70 to 70,000. Click on 
OK and then Solve. Stagnation is reached after 9.59 seconds on my laptop with the final results shown 
in Fig. 2.13 with SSR = 0.0842. 

Now we have the global minimum for the model equation of the boiling points we need to manually 
vary each constant in turn and look for the changes in SSR of a factor of 10. Such an approach quickly 
gives us estimates of the precisions of the arbitrary constants, A is precise to ±0.1; B to ±0.01; and C to 
±0.1. The NLLSQ fit of the model equation for the six n-alkanes is shown below and Fig. 2.14 shows the 
data and our fitted NLLSQ equation. 


BP = (724.5± 0.1K)log(« + (4.07± 0.01))- (384.7± 0.1 K) 



Figure 2.14: boiling points of n-alkanes. 
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The boiling point of n-nonane is interpolated from our NLLSQ equation as, 

BP= (724.5± 0.1 K)log(9 +(4.0710.01))- (384.71 0.1 K) 
BP= (424.0 t 424.5 i 423.6)K 


BP = 424.01 0.5 K 


Our NLLSQ equation may be used to interpolate the missing boiling points for n-7 and n = 9 with 
a confidence of ±0.5 K. However, extrapolating outside the range of n values for other alkanes is much 
less certain and should be treated with extreme caution. 


Return to Question 1 (see page 38) 


2.3.2 Solution 2: Heat Capacity of Solid Lead 
The Solver spreadsheet for the Shomate equation is 



A 

B 

C 

D 

E 

F 

G 

1 

77 K 

t = 771000 

C(p)/J K 1 mol- 1 

C(p calc) 

SR 



2 

300 

0.30 

26.7 

26.71 

0.00019 

A = 

26.839 

3 

350 

0.35 

27.3 

27.26 

0.00190 

B = 

2.028 

4 

400 

0.40 

27.7 

27.72 

0.00039 

C = 

4.000 

5 

450 

0.45 

28.1 

28.15 

0.00229 

D = 

0.916 

6 

500 

0.50 

28.6 

28.56 

0.00125 

E = 

-0.101 

7 

550 

0.55 

29.0 

28.98 

0.00026 



8 

600 

0.60 

29.4 

29.41 

0.00019 



9 








10 




SSR = 

0.00647 




Figure 2.15: Shomate equation for solid Pb C p NLLSQ fitting. 


Column D has the calculated heat capacity with D2 having, 


= $G$2 + ($G$3*B2) + ($G$4*(B2 A 2)) + ($G$5*(B2 A 3)) + $G$6/(B2 A 2)) 


which refers to the temperature in kelvin/1000 and the location of the five arbitrary constants to A to 
E. Cell E2 takes the residual and squares it, i.e. the E2 is =(C2-D2) A 2. The next question is the starting 
estimates for the constants, typing in unity into each of them gives SSR = 3037. Looking at the input data 
and the model equation then the heat capacities are in the 20s, so make A = 20 and typing in some values 
for B we find an initial value for the first temperature term B = 5. This gives SSR = 63.4 when the other 
arbitrary constants are unity, this is a good enough starting point for iterating to the global minimum. 
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With the cursor on the SSR click on Tools, Solver, and in the Solver window, minimize and select 
cells G2 to G6 to vary. Click on Options and make sure the solver engine is DEPS; edit the number 
of learning cycles (iterations) from 2,000 to 2,000,000 and increase the stagnation number from 70 to 
70,000. It took 15.7 seconds and 1454 iterations to find the fit with SSR = 0.00647 which is a good fit, 
see Fig. 2.16. Finding a measure of the precision of the constants is carried out by changing the digits 
(and then undoing the change) systematically through the each constant until we find a change of about 
an order of magnitude in the SSR. This gives us a final Shomate equation which the final figure shown 
is the last significant digit. For clarity I have not included units in the equation. Fig. 2.16 shows a graph 
of the experimental data with our fitted line for the above Shomate equation. 

C p = (26.8± 0.1 )+(2.0± 0.1)/ + (4.0± 0.1)/" + (0.9± 0.5)/' - ( al0± ; 0 01 ) 



Empowering People. 
Improving Business. 


Norwegian Business School is one of Europe's 
largest business schools welcoming more than 20,000 
students. Our programmes provide a stimulating 
and multi-cultural learning environment with an 
international outlook ultimately providing students 
with professional skills to meet the increasing needs 
of businesses. 


B! offers four different two-yea i; full-time Master of 
Science (MSc) programmes that are taught entirely in 
English and have been designed to provide professional 
skills to meet the increasing need of businesses.The 
MSc programmes provide a stimulating and multi¬ 
cultural learning environment to give you the best 
platform to launch into your career 

* MSc in Business 


* MSc in Financial Economics 


* MSc in Strategic Marketing Management 


NORWEGIAN 
BUSINESS SCHOOL 


EFMD 

EQUIS 


*ffi 


* MSc in Leadership and Organisational Psychology 

www.bi.edu/master 



Download free eBooks at bookboon.com 











Week 2: Chemical Data: Non-Linear 

Advanced Maths for Chemists Least Squares Curve Fitting 



Return to Question 2 (see page 39) 


2.3.3 Solution 3: Ionization of He 

Fig. 2.17 is the spreadsheet for the NLLSQ fitting of the ionization data for helium for our model equation. 



A 

B 

C 

D 

E 

F 

1 

E/eV 

S/A 2 

S(calc) 

SR 



2 

150 

0.43 

0.420136 

0.000097 

A = 

0.446103 

3 

175 

0.40 

0.407121 

0.000051 

B = 

426.554 

4 

200 

0.39 

0.393236 

0.000010 



5 

250 

0.37 

0.365113 

0.000024 



6 

300 

0.34 

0.338473 

0.000002 



7 

350 

0.31 

0.314232 

0.000018 



8 

400 

0.28 

0.292532 

0.000157 



9 

450 

0.27 

0.273214 

0.000010 



10 

500 

0.26 

0.256024 

0.000016 



11 

550 

0.24 

0.240695 

0.000000 



12 

600 

0.23 

0.226982 

0.000009 



13 

650 

0.21 

0.214665 

0.000022 



14 

700 

0.21 

0.203558 

0.000041 



15 

750 

0.19 

0.193502 

0.000012 



16 

800 

0.19 

0.184362 

0.000032 



17 

850 

0.17 

0.176022 

0.000036 



18 

900 

0.17 

0.168387 

0.000003 



19 

950 

0.17 

0.161372 

0.000074 



20 

1000 

0.16 

0.154907 

0.000026 



21 







22 



SSR = 

0.000642 




Figure 2.17: ionization cross-sections for helium electroionization. 


Column C is the calculated cross section S e.g. cell C2 is, 

= ($F$2) - ($F$2*EXP(- $F$3/A2)) 
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and the SR in D2 is =(B2-C2) A 2. The formulae from C2 and D2 are copied down to C20 and D20 with 
the appropriate variable for column A. The SSR in cell D22 is =SUM(D2:D20). 

The next important task is the initial guess for A and B. Entering the default A = 1 and B = 1 gives the 
sum of squares equals to 1.42. The first value for S is around 0.4 so A should be around 0.4. Trial and 
error by typing in various values for the constants with A = 0.4 and typing in B as 100, 200, 300, 400, 
500 and 600 gives shows that B = 500 has SSR = 0.004 which is the smallest value. So the initial estimates 
are A = 0.4 and B = 500. Selecting DEPS and editing cycle limit to 2,000,000 and stagnation limit 70,000. 
Solver reaches stagnation after 17 seconds with 1096 iterations and SSR = 6.42xl0“ 4 . 

Altering the digits of each constant looking for changes of an order of magnitude in SSR gives the 
following model equation. The NLLSQ equation is shown below with the last significant figure shown 
A = 0.45 ± 0.01 A 2 and with B = 430 ± 20 eV. Fig. 2.18 shows a graph of the experimental data with our 
fitted equation. 


S(E)= (0.45A 2 )- (0.45A 2 )exp-(-|^) 



Return to Question 1 (see page 39) 


2.3.4 Solution 4: Enzyme Reaction Kinetics 

In Fig. 2.19 the enzyme catalyzed data are plotted, along with a curve which is v = 0.21 [S]/(18 + [S]) with 
the estimated constants obtained by “eye” but nevertheless the curve is a reasonable initial fit to the data. 
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Thus the initial guesses are v max = 0.21 mol L 1 s 1 and K M = 18 mol L 1 and the NLLSQ solved spreadsheet 
is shown in Fig. 2.20. 



49 

Download free eBooks at bookboon.com 










Advanced Maths for Chemists 


Week 2: Chemical Data: Non-Linear 
Least Squares Curve Fitting 



A 

B 

C 

D 

E 

F 

1 

[ATP] 

V 

v calc 

SR 



2 

7.5 

0.06 

0.068014 

0.000064 

vmax = 

0.219857 

3 

12.5 

0.10 

0.093976 

0.000036 

£M = 

16.743864 

4 

25 

0.13 

0.131670 

0.000003 



5 

32.5 

0.15 

0.145102 

0.000024 



6 

50 

0.16 

0.164702 

0.000022 



7 

75 

0.18 

0.179732 

0.000000 



8 

100 

0.19 

0.188324 

0.000003 



9 

125 

0.19 

0.193886 

0.000015 



10 

150 

0.20 

0.197780 

0.000005 



11 

175 

0.21 

0.200658 

0.000087 



12 

200 

0.20 

0.202873 

0.000008 



13 

225 

0.21 

0.204629 

0.000029 



14 

250 

0.20 

0.206057 

0.000037 



15 

275 

0.21 

0.207239 

0.000008 



16 

300 

0.21 

0.208235 

0.000003 



17 

325 

0.21 

0.209085 

0.000001 



18 

350 

0.20 

0.209820 

0.000096 



19 







20 



SSR = 

0.000441 




Figure 2.20: Michaelis-Menten equation NLLSQ fit of APT-myosin data. 


The worth of the approximate graph of the data is shown when the v max and K u are made unity the SSR 
is 10.95. However, when are estimates from the Fig. 2.19 are typed in (v max = 0.21 and K u = 18) then 
SSR = 0.0021. Using the DEPS algorithm and editing the cycle limit to 2,000,000 and stagnation limit to 
70,000 Solver reaches stagnation after 14.7 seconds with 1088 iterations and SSR = 0.000441. 

Manually altering the maximum rate of reaction v max and Michaelis constant K M each one by one digit 
at a time (and undoing the changes) we look for the last significant figure (defined as altering the sum 
of squares of the residuals by about an order of magnitude). The two constants are then found to be, 


v max = (0.22± 0.01 mol L~ 1 s -1 ) and= (l7±4molL _1 ) 


and the Michaelis-Menten Equation is below with the last significant figures shown, 


(0.22 mol L~ l s~ l )[S mol L" 1 ] 
(l7molL 1 ) + [S molL" 1 ] 


The NLLSQ fitted plot is in Fig. 2.21. 
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Clearly with the linear and non-linear methods of curve fitting we have covered in Weeks 1 and 2 you 
can now tackle any model equation and experimental data in the future, including when you are doing 
research. 


Return to Question 2 (see page 41) 



Brain power 
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3 Week 3: Chemistry and Complex 
Numbers 

3.1 Why use Imaginary Numbers and Complex Numbers? 

I am giving a longer than normal introduction to imaginary and complex numbers because as a student 
I couldn’t see why some lecturers and professors wanted to use complex numbers instead of tangents, 
sines and cosines which I knew from school and also because the explanations in Maths books are not 
always very helpful. I want to show you why they are so useful, they cut down the use of a calculator 
and reduce the use of trig relationships. I want to start by comparing imaginary numbers with their 
close neighbours, the negative numbers. Negative numbers give you a “mirror image” of the positive 
number on the other side of zero on the number line, thus positive and negative numbers are all one¬ 
dimensional. Multiplying x by -1 “flips” you along the number line to -x. 

real dimension 


Figure 3.1: the positive and negative real number line. 

For example in words, “I owe the bank £1000 and after doing some evening work I earn £50 and then 
pay the bank £50, then I now owe the bank £950”. A better, and more elegant way of handling the 
situation is to use positive and negative numbers and Maths which keep track of the direction along the 
number line automatically my money = -£1000 +£50 = -£950. You shouldn’t be too smug about the 
use of negative numbers (don’t say “I know all this, so why bother”, well be patient). Negative numbers 
were not understood and were considered absurd by many great Mathematicians until the 1700s then 
they only came into normal usage driven partly by their short elegant way of summarizing things from 
book keeping to Chemical Analysis. 

Imaginary numbers have a similar history to negative numbers, in that people didn’t or still don’t think 
they are “real” (pun intended). We can solve for x for equations such as x 2 = 25 easily, but what does 
the solution for x 2 = -25 mean? Although they may appear strange to us, they are no more strange than 
zero, negative, fractional and irrational numbers (which can’t be written as a simple fraction) appeared to 
scientists before us when they were first introduced. Let us assume the imaginary number i exists, where 

i = (“ 1) 12 
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As i is just a number (like e and tt) it is written in upright roman text. If i exists then Maths becomes 
easier and new relationships emerge that are easily described. I dislike the term “imaginary number” 
it was introduced when i was discovered around 1545 by the first group of scientists who thought they 
were nonsense and absurd. The number i is just as normal as other numbers but the name “imaginary” 
stuck for purely historical reasons and the terms “real” and “imaginary” are just labels. 

3.1.1 Visual Understanding of Negative and Complex Numbers 

The equation x 2 = 25 really means “what operation on v, when applied twice , turns 1 to 25” i.e. 

lXx 2 = 25 

The two answers are x = 5 and x = -5, i.e. you can “scale 1 by 5” or “scale 1 by 5 and then flip” and do 
either of these operations twice over, where “flipping” or taking the opposite is multiplying by a negative. 
We are still dealing with the one-dimensional real number line. 

What about the imaginary number x 2 = -1? 


IXx 2 =-l 

“What operation on v, when you apply twice will turn 1 into -1”? It cannot be multiplying by a positive 
twice because the result remains positive. It cannot be multiplying by a negative twice as the first negative 
multiplication flips v to negative but the second one flips it again to give a positive number. So the 
operation we are looking for cannot be just on the real number line of positive and negative numbers 
(labelled “Re”). 

We need to add a second dimension (labelled “Im”) at right angles to the real number line for imaginary 
numbers which crosses at the origin (0,0), Fig. 3.2. We are now dealing with two-dimensional number 
space called the complex plane and Fig. 3.2 is an example of an Argand diagram named after Jean-Robert 
Argand (1813) who first popularized the idea of visualizing imaginary numbers in terms of a graph 
although people had tried to understand complex numbers from around 1545 onwards. 
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So returning to “what operation on when you apply twice will turn 1 into -1” the operation we are 
looking for is anticlockwise rotation of 1 by 90° takes % onto the imaginary line at i and then another 
anticlockwise rotation of i by 90° takes i to -1. Alternatively, we could rotate twice in the clockwise 
direction to turn 1 into -i and another clockwise rotation turns -i into -1. So there are two square roots 
of -1, they are i and -i. 


x = 


<-ir=±i 
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3.1.2 Patterns of Numbers 

For negative numbers, if we repetitively multiply xby -1 we get the pattern below, 

x, — x, x, - x, x, •• 


because -1 does not change the size of a number only its sign we just flip back and forth along the real 
axis. So negative numbers model things which flips back and forward (oscillate) between two situations. 

For imaginary numbers, if we repetitively multiply by i we get the pattern below, 

x, ix, - x, - ix, x, - 

this is easy to visualize as rotations on an Argand diagram, Fig. 3.3. We have completed a 360° cycle 
after every fourth rotation. 



So although negative numbers model oscillations along the real axis, imaginary numbers model rotations 
in two-dimensions i.e. anything with a cyclic or a rational motion. So wave motions as in wavefunctions 
are ideal candidates as well as other cyclic processes that Chemists have to deal with. But why not use 
tan, cos, sin instead?. The answer is it is easier, less calculator work, less pure memory work of trig 
relationships, we cover this aspect now. 

3.1.3 Complex Numbers 

A number may be both “real” and “imaginary” and we call it a complex number. We may rotate by any 
angle not just 90°. So if we rotate a number x by any arbitrary angle 9 other than 90°, the result would 
have both real and imaginary parts and be a complex number which is written as 


z = a + ib 
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Where a and b are both real numbers, but i b is the imaginary part of z and a is the real part. The 
components a and b may be either positive or negative, i.e. the complex number may be in any of the 
four quadrants of the Argand diagram. At first the idea that a single number has two parts is a bit strange 
but you have been doing something similar for years. What does the decimal 10.7 mean? It consists of 
two parts, an integer and a fraction and so 10.7 = 10 + 7/10, i.e. one number consisting of two parts. The 
angle 0 is called the argument, arg(z) and is the angle of the complex number measured from the real 
axis, with anticlockwise rotation defined as a positive angle. The word “complex” just means “a mixture” 
when used to describe numbers. 

3.1.4 Size of a Complex Number 

The size or magnitude of a negative number, -x, is its absolute distance from zero (called the modulus) 
and written as |x|. 


\x\ = V—x 2 

The size of a complex number, its modulus |z|, is the absolute distance from zero which we can obtain 
using Pythagoras’s theorem in Fig. 3.4 and is the absolute length of the red line. 



3.1.5 Forms of a Complex Number 

z = a + ib Standard Form 

The above is the standard form. From Fig. 3.4 we have a = \z\ cos(0) and b = \z\ sin(0) which gives the 
polar form below, where 0 may be in degrees or radians. 


z=\z\ (cos 0 + i sin 0) Polar form 
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From the definition of an exponential as a power series and then replacing v by id we obtain, 


exp(i0 ) = 1 + (i0) 


2 
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from the series expansions for sin and cos we have Euler’s formulae for positive and negative exponents 
and which then gives the exponential form of a complex number. Where 9 must be in radians so that 
W of the exponential is dimensionless. 


exp(i0)= cos0 +isin0 exp(-i0)= cos0 - isin0 Euler's Formulae 


z=\z\ exp (i0) Exponential form 


From Eulers formula (1748) we can get the log version of a complex number. But the important thing is 
there is very little memory work (probably best to remember Euler’s formulae) if we use the geometric 
approach of the Argand diagram which allows easy conversion between the various forms. 

iG = ln( cos0 + isin0) Logarithmic Form 
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3.1.6 Complex Addition 

Addition of real numbers may be thought of as “sliding” along the real number axis. Addition with 
complex numbers is similar, but we “slide” in two dimensions simultaneously but separately (we add the 
real and imaginary components separately). For z x and z 2 shown below then the sum z 3 is, 

z l = (2 + i 4) z 2 = (—3 + i 2) 
z 3 = (2 + i4)+ (- 3 + i2)= (- 1 + i6) 

This is easier to understand from the Argand diagram (Fig. 3.5) which shows the “sliding” in the two 
dimensions with z x , z 2 , z 3 coloured (black, red, and dashed green for clarity). 



Figure 3.5: complex addition. 


3.1.7 Complex Subtraction 

The process of subtraction is the “opposite” or reverse process to addition. As before, for z and z 2 , then 
the difference z 3 is equal to, 


z x = (2 + i 4) z 2 = (-3 + i 2) 

z 3 = (2+i4)- (- 3 + i2)= (2 + i4)+ (3 - i2)= (5 + i2) 

In the Argand diagram, we draw z x and the negative of z 2 (where -z 2 has the same magnitude as z 2 but 
the real and imaginary parts change sign and so the angle is at 180° to the original z 2 ). We then add 
z x and -z 2 following the normal rules for adding the two components, see Fig. 3.6. 
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Re 


Figure 3.6: complex subtraction. 


3.1.8 Complex Multiplication 

Multiplication and division are the “opposites” or reverse operations of one another. So starting with 
multiplication it is simply multiplying the two complex numbers using the normal rules of algebra for 
multiplying out a pair of brackets, as below. 

z, = (2 + i4) z 2 = (—3 + i 2) 

z 3 = (2 + i 4 X - 3 + i 2) = - 6 + i 4 - i 12 +i 2 8 


z = — 14 — i 8 


9 , = tan -1 (^1 = 29.7° 


The complex multiplication above is actually carrying out the following two operations, 


Complex Muliplication: muliply the sizes and add the angles 


kl =kl*kl 


0 , = *. + *2 


Also notice that we really multiply the sizes when we multiply real numbers, 3x4= |3|x|4| = 12 and also 
with negative numbers, -3x4 = -(3x4) = —(|3|x|4|) = -12. As before, for z x and z 2 as given preciously, 

IzJ = V 2 2 +4 2 = V 20 |z 2 |= V— 3 2 + 2 2 = Vl3 


z 3 | = V20xVl3 = 16.1 


This agrees with the size from z 3 = V(14 2 + 8 2 ) =16.1245 for our straight multiplication. The angle is the 
sum of the angles, 
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In the Argand diagram, Fig. 3.7, z 1 = (2 + i2) is in the first quadrant; z 2 = (-3 + i2) is in the second 
quadrant and the angle d 2 is negative (clockwise) from the real axis. The resultant z 3 = (-14 - i8), not 
shown in full for clarity, is in the third quadrant and the angle d 3 is positive (anticlockwise) from the 
real axis. 



Figure 3.7: complex multiplication. 


3.1.9 Complex Division 

Multiplication and division are the “opposites” or reverse operations of one another. 
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How do we subtract the angle? We use the complex conjugate z* which has the same real component 
but the imaginary component is the mirror image (the real axis acting as the “mirror”) as in Fig. 3.8. 



So multiplying by the complex conjugate is the same as subtracting the angle. So far we are at the stage of, 


z 2 

which has subtracted the angle all right, but at the same time it has altered the size of the numerator! 
So in order the correct this size change in the numerator we must alter the size of the denominator by 
the same amount, 



For the division of, 


z l = (2 + i 4) z 2 = (-3 + i 2) 

= ( 2 + i 4) (— 3 - i2) _ -6- i4- i!2- i 2 8 _ 2 - i 16 = 2 .16 
3 (-3+i2)(-3- 12) 9 + i 6 -i 6 -i 2 4 9 + 4 13 1 13 


z 3 = 0.1538 - i 1.2308 


Which gives the resultant angle (argument) of, 


= tan 


-1.2308 
i 0.1538 


-82.9° 


z 3 has a positive real part and a negative imaginary part so it is in the fourth quadrant and the angle is 
clockwise (negative) from the real axis. This is clearer from the Argand diagram Fig. 3.9, z 3 is not shown 
in full for clarity. 
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Figure 3.9: complex division. 


3.1.10 de Moivre's Formula 

Multiplying the polar form of the complex number, 

z = \z\ (cos 9 + i sin 9) 

by itself n times then we multiply the size |z| by itself n times, giving the final size \z\ n . We add the angle 
9 to 0, n times, giving the final angle n9. The resultant polar form is called de Moivres formula (also 
known as de Moivres theorem), 


z n = \z\ n (cos n9 + i sin n9) 


where n may be positive, negative, integer, or fractional, although if n is fractional then the result is 
multi-valued. I will give the two most common examples of the use of de Moivres theorem, the square 
root and the cube roots of of a complex number, 

z = 4 (cos 60° + i sin 60°) 


As z is identical for 60° and 60° plus one revolution i.e. for 60° + 360° = 420°, we have, 




1/2 .i/ 2 # 60° . . 60° 

z = 4 I cos —-—y l sin—— 


) 


(‘ 


, 1/2 . 1 / 2 # 420°. . 420° 

and z =4 I cos—— + 1 sm—— 


) 


z 1/2 = 2(cos30°+ i sin30°) and z m = 2(cos210° + i sin210°) 


There are two square roots of a complex number separated by 180°, i.e. 360°/2, as in Fig. 3.10. All you need 
to do is find the first root (the principal root) and then add another root at 180° to the first square root. 
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Figure 3.10: square root of a complex number. 


The second example of the n being fractional is the cube root of, 

z = 27(cos 120° + i sin 120°) 


As z is identical for 120° and 120° plus either one or two revolutions i.e. (120° + 360° = 480°) and (120° 
+ 360° + 360° = 840°) we have, 




1/3 _ ryrj 1/3 i 120 • • 120 

z =27 (cos-+ l s in- 


) l/3 _ 

z - 27 Icos 


480' 


. 480° 
— + i sin—— 
3 3 


) 


1/3 


1/3 


z =27 (cos 




840° 


840° 


„ + i sin „ 
3 3 


) 
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z 1/3 = 3(cos40° + i sin40°) z 1/3 = 3(cos 160° + i sin 160°) z 1/3 = 3(cos280° + i sin280°) 


There are three cube roots of a complex number each separated by 120°, i.e. 36073, as in Fig. 3.11. All 
you need to do is find the first root and then add 120° twice for the other two cube roots. 



3.2 Week 3 Tutorial Questions 

3.2.1 Question 1: Quantum Chemistry 

In the quantum mechanics of electrons, atoms, molecules, or ions; chemists make use of wavefunctions 
(if/). The reason is that the value of the wavefunction often varies in a cyclic fashion and so can very 
naturally be expressed as a complex number. If the wavefunction f has the value 

ijj - A + iB 

the probability of the system having that value of f is related to the product of the complex conjugate 
i//* and the wavefunction itself \)/* f. Find whether the probability of the system is a real number? 

Tump to Solution 1 (see page 66) 


3.2.2 Question 2: Atomic Orbitals 

The wavefunctions for the 2p x and 2 p atomic orbitals for the H-atom can be written as, 


V 2p = A eX P 


— r sin0 exp(i^) y/ 2 = A exp — — ^ — r sin# exp(—i cp) 

2 a o] y \ 2 a o) 


Notice that the 2 p x and 2p y atomic orbitals only differ in the sign of the complex number term for the 
phi exponential. Fig. 3.12 shows the polar coordinates (r, 0, (p) where r is the radial distance from the 
origin at the centre of a sphere; 9 (theta) varies from 0 at the “north pole” to n radians at the “south 
pole”; and cp (phi) varies around the “equator” from 0 to 2n radians. 
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Figure 3.12: polar coordinate system. 

Rewrite the sum of the two wavefunctions for 2p x and 2p^ atomic orbitals for the H-atom as a real number. 

Tump to Solution 2 (see page 67) 


3.2.3 Question 3: Chemical Bond Vibrations 

Chemists use infra red spectroscopy (IR) which relies on the vibrations of chemical bonds to identify 
particular “functional groups” in an organic or inorganic compound and as a method of identifying 
unknown molecules. Fig. 3.13 shows the bond stretching and shortening in a diatomic molecules bond 
length with time from a fully compressed distance, through the equilibrium (lowest energy) distance, 
to the fully extended bond length and then all the way back again. 

0—0 0rO 

fully compressed equilibrium fully extended 

bond distance bond distance bond distance 

Figure 3.13: simple harmonic vibration of a chemical bond. 

The theory behind this important spectroscopic technique makes use of a simple mathematical model to 
describe the vibration of a chemical bond which is called simple harmonic motion (SHM). The length 
of a chemical bond undergoing SHM can be described by the following equation. 

r[t) — r e + A cos [kcot] + B sin [hoot] 
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Where t is time; co is the angular frequency in radians per second of the vibration with co = 2nv where 
v is the ordinary frequency in hertz (s _1 ); A and B are arbitrary constants; r e is the equilibrium bond 
distance; and k is the force constant or “stiffness” of the chemical bond. 

The periodic motion of the atoms in a vibration are more naturally modelled using complex numbers 
rather than trig functions as the bond length varies in a cyclic fashion with time. Write down an alternative 
expression for this SHM equation in term of a complex exp (i kcot) and its complex conjugate. 

Tump to Solution 3 (see page 68) 


3.3 Week 3 Tutorial Solutions 

3.3.1 Solution 1: Quantum Chemistry 

Multiplying the wavefunction by its complex conjugate is easily done with basic complex number maths. 

\f/= A+iB 

y/*y/= (A - i B)(A + i B)= A 2 + iAB - i BA - i 2 B 2 


y/* y/ = A 2 + B 2 
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The product of y* f is always real and always positive (> 0), i.e. a wavefunction probability at a given point 
in space will be a real positive number, as it should be if the concept of probability is to make sense. Fig. 
3.14 shows the wavefunction (red) for the n-2 quantum level of the particle in a 1-dimensional box is 
both a negative and positive complex number where the green line is the zero level for both f 2 and y* 
y r The probability of the particle being at a certain position (blue) in the box is a positive real number. 



Return to Question 1 (see page 64) 


3.3.2 Solution 2: Atomic Orbitals 

The wavefunctions for the 2 p x and 2 p atomic orbitals for the H-atom can be written as, 


v 2p = A ex P 


2 a n 


r sin # exp (i q>) yj 2p — A exp 


2 a n 


r sin# exp — icp) 


To rewrite these as purely real numbers we make use of Euler s formulae written in terms of i<p, 
exp(i^) = cos cp + i sinyp and exp(— icp) — cos cp — i sin^9 

Adding our two wavefunctions, 


¥ 


2 p 


+ Vip = A exp 


r 
2 a 


r sin# exp i</9 + exp — \cp 


°/ 


Substituting using Eulers formulae gives the sum of the 2p x and 2atomic orbitals for the hydrogen 
atom which is a real function. 


yj 2 p +y/ 2p = ^4 exp I— -— I r sin# (cos<p + i sing9 + cosq 9 — i sin(/9j 


^ 2 p + ^ 2 P = 2 A exp 


2a °, 


r sin # cos op 
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In quantum chemistry, it is often more convenient to use complex-functions, particularly when one needs 
to apply quantum mechanical operators to a wavefunction to generate a new wavefunction. 

Return to Question 2 (see page 64) 


3.3.3 Solution 3: Chemical Bond Vibrations 

r( t)= r e + Acos(kut) + Bsm(kU)t) 

exp(i0) = cos0 + i sinS and exp(- id ) = cos0 - i sinS 

Euler’s formulae for positive and negative exponentials may be added together and after rearrangement 
gives cos(0) in terms of complex numbers. Subtraction and rearrangement gives sin(0) in terms of 
complex numbers. Note that we only need to remember Euler’s formulae we do not need to remember 
the trig equations, certainly I don’t remember them. 

cos0 = (exp(i@ ) + exp(- i@)) and sin0 = 7 y(exp(i@)- exp(-i@)) 


We can now substitute for cos ( kcot ) and sin ( kcot) in our equation for simple harmonic motion where 
we have let 9 = kcot. 

A B 

r(t)= r Q + — (exp(i^cu^) + exp(- ikoot)) + — (exp(i^6UZ L ) - exp(- iku) t)) 


Expanding out the brackets, 


A A B B 

r(t)= r Q + — exp(i^6U t) + — exp(- i kcot) + — exp(i£cu t) - — exp(- ikoot) 

A* A* 1 A* 1 


In the last two terms we are diving by i which is the same as multiplying by -i which is shown below as 
an aside and comes from multiplying and dividing by i, 


1 

i 



i 


i 

i 



- l 


A A iB iB 

r(t)= r e + — exp(i^cu^) + — exp(- i kcot)- —Qxp(iku)t)+—Qxp(-iku)t) 


Collecting the terms with the same exponents gives, 


/x A - iB .. 7 x A + iB . . 7 x 

r(t)= r e +— --exp(i^6U^) + —--exp(-i^6U^) 


Return to Question 3 (see page 65) 
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4 Week 4: Chemistry and Complex 

Numbers 

4.1 Week 4 Tutorial Questions 

4.1.1 Question 1: Particle In A One-Dimensional Box 

Quantum effects are dominant in many areas of Chemistry and Biology A simple system that shows 
how quantum effects arise is the particle in a box. The particle may be an electron, atom, ion or molecule 
dependent upon the system being studied. The first three wavefunctions for n = 1,2, and 3 are shown in 
Fig. 4.1 which has a vertical axis not to scale. The wavefuction increases and decreases with the distance 
along the box and thus is best modelled using complex numbers. 

The wavefunctions f for the particle in a one dimensional box are given by the following complex-function 

= j~2 (exp(i^) - exp(-iwx)) 
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Where L is the length of the box, x is the position of the particle within the box and n is a constant called 
the translational quantum number n- 1,2,3,4, Simplify this expression to obtain a real valued function. 



^3 

^2 


Figure 4.1 : lower energy wavefunctions for a particle in a box in red, with zero values in green. 


Tump to Solution 1 (see page 73) 


4.1.2 Question 2: X-Ray Crystallography Structure Factors 

X-ray diffraction is a standard method of chemical analysis to determine the structure of crystalline 
solids. The X-rays photons are scattered from the electrons in the atoms and interfere to give bright spots, 
i.e. the crystal planes of atoms act as diffraction gratings. The crystal planes of atoms are described by 
Miller indices (h k l ) from the ratios of the planes coordinates (x, y, z) to the lengths of the axes of the 
unit cell a, b , c (which are not necessarily at right angles). We then take the reciprocal of these ratios; 
and if necessary multiply these reciprocals by the smallest number to converted them to integers. In Fig. 

4.2 the crystal plane are shown in blue and these Miller indices become really useful for the planes that 
are parallel to the an axis. So for the top diagrams the first plane which is not at the origin is a la, and 
it is parallel to the b and c axes so these intercepts are both at la, oo b, oo c which gives the ratios and the 
planes Miller indices using the fact that l/°o = 0. 

f — — {hkl)= (100) 

1 OO oo 

In the middle diagram the first plane not at the origin has intercepts at Via, Vib , and oo c, giving the ratio 
and Miller indices, 


11± (hkl)=( 220) 

/2 /2 oo 

In the bottom diagram the first plane not at the origin has intercepts at a, b , c which gives the ratios 
and Miller indices. 


yjy (**/)=( HI) 
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(1 0 0) planes 


(2 2 0) planes 
c 

(111) planes 


Figure 4.2: a crystal with two unit cells and three types of crystal planes. 

The intensity of the scattered beam of X-rays from the hkl plane of a crystal is proportional to FT where 
F, the structure factor for a particular crystal is determined by the positions of the atoms in the various 
crystal planes and F is its complex conjugate. 




F {hkl)- ^ f j exp(i2 TT {hx j + kyj + lz j)) 

j 


The summation runs over j the number of atoms in the unit cell, with each atom having a scattering 
factor^ and fractional coordinates (x., y. y z.) within the unit cell. 

Metallic sodium crystals have a body-centred cubic structure with two atoms per unit cell located at 
(0,0,0) and (V 2 F/ 2 F/ 2 ) respectively. Firstly, use the above equation to find F the structure factor for metallic 
sodium, and secondly, use Euler s formula to determine its real and imaginary parts. Thirdly, show that 
X-ray diffractions “reflections” occur (i.e. F(hkl) ^ 0) only if the quantity h+k+l is even, but not if it is odd. 

Tump to Solution 2 (see page 74) 

4.1.3 Question 3: Dielectric Constant of a Solution of Polar Molecules 

An alternating electric field of frequency co has an electric field, E = E Q exp(i cor), which varies with time 
r (“tau”). Applying an alternating field to a dilute solution of polar molecules allows you to measure the 
effect of the solvent with the polar molecules in the solution. This solvent effect is called the dielectric 
constant e of the solution at the frequency co of the applied field. For example, water which is very polar 
has a dielectric constant of 80.10 at 20°C but n-hexane which is very non-polar has a dielectric constant of 
1.89 at 20°C. This information is of great value in deciding which solvent(s) to use in designing separation 
methods, for sample preparation and in the chromatography techniques used in analytical chemistry. 

The dielectric constant of a solution is, 


e - 


2 

n = 


jVoPA ' 2 

3M e 0 k B T (1 + i out) 
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The dielectric constant e is a complex number and is conventionally written as, 

£=£'+{£" 

The parameters n , N {) , p , M, a 0> p, k B , T are all constants at a constant temperature T. Experimentally we 
measure the real part of the dielectric constant e so rewrite the above equation to determine the real 
part e! of this complex dielectric constant. 


Tump to Solution 3 (see page 75) 

4.1.4 Question 4: Polymer Electrolytes and Impedance Spectroscopy 

New batteries for car, laptops, tablet computers and mobiles are developed by Chemists. These batteries 
are now made from polymer electrolytes as these can offer significant performance improvements. 
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One tool that is often used to characterise the behaviour of these polymer electrolytes is the technique 
of impedance spectroscopy, where the electrical impedance Z of a sample is measured as a function of 
an applied alternating current of angular frequency co. The impedance Z is the opposition to the flow 
of an alternating current and is measured in ohms. Angular frequency co = 2n/T where T is the period 
of rotation in seconds, so co is in radians/second and is equals 2nv where v is the ordinary frequency in 
hertz or s 1 . It is found experimentally that some materials exhibit an impedance response that is given 
by the following equation 


Z(CJ) = 


R 


o 


1 + i COT 


Where R is the direct current (DC) resistance of the polymer electrolyte and is a constant and r (tau) 
is the time. It is found that impedance is a complex variable, 

Z( UJ)= x(GU)+iy(GU) 


Firstly, find the expressions for the real and the imaginary parts of the impedance equation i.e. find 
expressions for both x(co) and y(co). Secondly, eliminate cor from these equations and graph Z(co) in the 
complex plane as a function of co and show it is a circle. Find the (x, y) coordinates of its centre and 
the radius of the circle. For a good polymer electrolyte for a new battery, we will needs a high ionic 
conductivity i.e. a low DC resistance, R Q . As a hint the equation for a circle with its centre at (a, b ) and 
radius r is, 

(x- af + (y- bf= r 2 

Tump to Solution 4 (see page 76) 

4.2 Week 4 Tutorial Solutions 

4.2.1 Solution 1: Particle In A One-Dimensional Box 


2 1 


V / \ x l = \j— — (exp^ircxj — exp^— inx, 


Substituting from Eulers formulae written in terms of (inx) we have, 

"l 1 


L i 2 


= y— — Icos^nxJ + isin(nxj — cos(nxJ + isin(nx, 


V\x) = \jj~ T i2sin(nx) 



This is a real wavefunction. 


Return to Question 1 (see page 69) 
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4.2.2 Solution 2: X-Ray Crystallography Structure Factors 

Firstly, from the definition of the structure factor, we have to substitute the values for each atom in the 
unit cell. 


F (hkl)= ^ fj exp(i2 TT (hxj + kyj + Iz y )) 


Atom 1 has a scattering factor/ Na and coordinates (0, 0, 0), 


/ Na exp(i2jc(/z0 + kO + /Oj) = / Na exp(i2jtOJ = / 


Na 


since exp(0) = 1 the first sodium atoms contribution to the structure factor F summation is just/ Na . 
Atom 2 has a scattering factor/ Na and coordinates (V 2 , V 2 , V 2 ) and makes a contribution to the summation of, 


/n a eX P 


i 2 jc 


| + | + | = / Na exp(ijt(/t + k + l) 


Adding the effects for the two atoms gives the structure factor F for metallic sodiums body-centred 
unit cell, 


F(hkl)= / Na + / Na exp(iTT(/* + k + /)) 


Secondly, using Eulers formula we will determine the real and imaginary parts of the structure factor F. 

exp(iTT)= cos(l80°)+ i sin(l80°)= - 1 + i(0)= - 1 

where we must use radians (tt radians = 180°). For the particular example of a body-centred cubic 
structure the real part is -1 and the imaginary part is zero, but this will not necessarily happen for other 
crystal structures! 

Thirdly, we have to investigate the behaviour of the function exp(fn:(/z + k + /)) for different values of 
h+k+l. As /z, k and l are all integer quantities their sum must also be an integer. Lets call this sum N = 
h + k +1 and we must now look at the behaviour of exp(iiVTt) for various values of integer N. So we have 
when N = 1, 3, 5, ... 


exp (i A tt) = cos (N 180°)+ i sin(A 180°)= - 1 + i(0)= - 1 

substituting this value for exp(iNTt) into F(h k l ) gives F(h k l) = 0 for N = odd. When N = 2, 4, 6 , ..., 

exp (i A tt) = cos (N 0°)+ i sin(A0°)= +1 + i(0)= +1 
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This means that the structure factor F is either of the two values given below dependent upon the Miller 
indices of the crystal plane, 

F(hkl) = 2 / Na for N = h + k + / = even 


F(hkl)= 0 for N= h + k +l = odd 


i.e. there is no X-ray scattering observed for h + k + / = odd crystal planes but only from h + k + l = 
even planes for the particular crystal structure of metallic sodium. 


Return to Question 2 (see page 70) 

4.2.3 Solution 3: Dielectric Constant of a Solution of Polar Molecules 

£- n= ^oRHl 

U 3M £ 0 k B T (\ + i cot) 


£' + £ " = 


n + 


N 0 PIJ 2 


3 M£ 0 k B T(l + hot ) 
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Consider the term 1/(1 + icor) as it is the only term containing i and the rest of the terms on the right hand 
side are just constants. The 1/(1 + i cor) term is division by a complex number and when we multiply it 
top and bottom by the complex conjugate of the denominator and simplify we find the real part of this is, 

1 _ 1 1 ~ ic or _ 1 ~ ic or 

1 + i COT l+ic OT 1-i COT 1+(j0 2 T 2 

1 ~ i CUT _ 1 _ • COT 

1 + CO T 1 + CO T 1 + CO T 

o 1 - 1 

1 + 1 U)T 1 + (jj 2 j 2 

So the real part of the original equation can be rewritten, see below, and it is this part of the dielectric 
constant of the battery polymer electrolyte that is measured experimentally 


£' = n + 


n 0 p/j 2 

3M £ 0 k B T (l + cu 2 t 2 ) 


Return to Question 3 (see page 71) 


4.2.4 Solution 4: Polymer Electrolytes and Impedance Spectroscopy 

This has a similar solution to the last question but with a twist in the plot! First evaluate the real and 
imaginary parts by multiplying top and bottom by the complex conjugate. 


Z< 0 l >) = 


R ,„ 


1 + icor 


Z(U)) = 


R 0 1 - 1U)T 

-^-x - 

1 + iU)T 1 “ 1U)T 


i . 2 2 

1 + 6 0 T 


Rr, 


i . 2 2 

1 +60 T 


i , 2 2 

1 + 6 0 T 


In impedance spectroscopy this is usually written as below which gives us the twist that y (to) is a function 
of x(co), 

Z (CO) = x(C0)+iy(C0) 

R o 

X (CO)— -~—7 

1 + 6u 2 r 2 


R 0 cot 

y(C0)= - 2-5 = “ (cot ) x(CO) 

1+6 0 T 


COT = - 


y(C0) 
X (CO) 


(COT 2 = 


, . 2 

y(C0) 

X (CO ) 2 
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Thus on substituting for (cor) 2 in the equation for the real part of the impedance, x(co ), and for clarity 
dropping the reminder that x and y are both functions of the angular frequency co we have, 



Now multiply both sides of this equation by the denominator (l+y 2 /x 2 )and simplifying further, 

2 

x + — = R 0 x 2 + y 2 = R 0 x x 2 - R 0 x + y 2 = 0 
x 


Add X AR 2 to each side so that all the terms can be written as squares (completing the squares method) 


^ » * 0 + 2 _ K 

x - R ox + — + y - — 



As suggested in the question, this is the equation for a circle with centre at (ViR Q , 0) and radius VzR 0 , with 
the circle lying in the x-y plane (the Re(Z)-Im(Z) plane) as in Fig. 4.3 which is drawn for R 0 = 10 ohms. 



Return to Question 4 (see page 72) 
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5 Week 5: Chemistry and Matrices 

5.1 Why use Matrices with Chemistry? 

One of the most obvious uses of matrices is to organize the data of several simultaneous equations in 
a compact form so it can be manipulated and easily solved. This is most commonly used in Analytical 
Chemistry A set of linear simultaneous equations consisting of numerical values for the coefficients A 
and the values b can be solved for the unknown variables x when it is written in matrix form Ax-b. As 
we need as many simultaneous equations as there are unknown variables we will have a square matrix 
for A. If we multiply by the inverse of A on both sides it gives A 1 Ax = A~ l b or lx = A~ l b where I is the 
identity matrix and so^ = A~ l b. So the problem is reduced to finding the inverse matrix. To find the 
inverse we have to do several things, firstly find the determinant of A, |A|, and secondly, the matrix of 
cofactors of A, C, and then the transpose of C (C 1 = adj A the adjunct of A). Here is a simple reminder 
of these steps for the set of simultaneous equations below we have the A matrix, 

5 x l + 2x 2 + 3 x 3 = 4 
4 x : + 1 x 2 + 6 x 3 = 3 
8x 1 + 0x 2 +7x 3 = - 2 
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A = 


5 

4 


1 6 
0 l) 


\a\ = +5(lX7 — 6 x 0 ) -2(4X7- 6X8) + 3 ( 4 x 0 - 1X8) 


\A\ = +5(7) — 2(—20) +3(—8) = +35 +40 -24 = 51 


M| = 51 


The matrix of cofactors, C, of A is 


+(1X7 

-6X0) 

-(4X7- 

6X8) 

+(4X0 

- 1X8)1 


+7 

+20 

-8 

-(2X7 

- 3x0) 

+(5X7- 

3X8) 

-(5X0 

-2X8) 

= 

-14 

+ 11 

+ 16 

+(2X6 

— 3 X 1) 

-(5X6 - 

3X4) 

+(5X1 

-2X4)1 


+ 9 

-18 

-3 


The transpose of C is CT is, 


C T 


7-14 9 

20 11 -18 
-8 16 -3 


Putting it all together, we have the inverse of A is now, 


A~ x 



20 11 -18 
\—8 16 -3 


The solutions of the simultaneous equations are as follows, 


x = A 1 b 


X ~\A\ 


Substituting in for the determinant of A, the transpose of cofactors (C the adjunct of A), and the values 
for b 


7 

-14 

9 

4 

20 

11 

-18 

3 

\- 8 

16 

_3 / 

H 


x 


51 


4(7) + 3(-14) + (—2)(9) 
4(20) + 3(ll) + (—2)( —18) 
4(-8) + 3(16) + (-2)(-3) / 


-0.627 


x = 


2.922 
0.431 t 


51 


-32 

149 

\ 22 / 
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Xj = —0.627 x 2 = 2.922 x 3 = 0.431 


Solving small numerical matrices by hand may be fine and it is something you need to do for exam 
questions ! But for large matrices I will show you later a simple spreadsheet approach. 

A second use of matrices in Chemistry is the efficient handling of the symmetry operations that apply to 
molecules. The shape of molecules is one of their most important properties and the shape is defined by 
the symmetry operations that they can undergo. A molecule has a particular shape or symmetry which 
may be transformed in three-dimensional space by rotations around its symmetry axes; reflections in 
its mirror planes; rotation-reflection (or improper rotation); inversions through a centre of symmetry; 
and the identity operation. 

A third area in Chemistry where matrices are used is the eigenvalues and eigenvectors of the determinant 
of square a matrix to find the allowed solutions to the quantum mechanics of the structures and reaction 
pathways of a molecule. This last area is covered in Week 6. You may not yet have covered some of these 
areas of Chemistry in your undergraduate course. Dont worry they will be introduced as you move 
through your Chemistry degree and you will already be able to do the Maths! 

5.2 Week 5 Tutorial Questions 

5.2.1 Question 1: Chemical Ore Analysis 

In the lab an analytical Chemist measures the fractional Pb, Fe, Cu and Mn content on a weight for 
weight basis for three ores, her results are shown in Fig. 5.1. 



Ore 1 

Ore 2 

Ore 3 

Pb 

0.1 

0.2 

0.3 

Fe 

0.2 

0.3 

0.3 

Cu 

0.6 

0.2 

0.2 

Mn 

0.1 

0.3 

0.2 


Figure 5.1: metal analysis of three ores. 


If the ores are mixed together so there are x ] kg of Ore 1, x 2 kg of Ore 2, and x 3 kg of Ore 3, the amount 
of each metal in the mixture can be computed as shown in Fig. 5.2. 


Amount of Pb = 

^Pb " 

O.lxj + 0.2x2 + 0.3x2 

Amount of Fe = 

^Fe = 

0.2xj + 0.3 X 2 + 0.3 x^ 

Amount of Cu 

o 

c 

II 

0.6xj + 0 . 2 x 2 + 0 . 2 x 2 

Amount of Mn = 

^Mn~ 

O.lxj + 0.3x2 + 0-2x2 


Figure 5.2: ore mixture. 
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Firstly, rewrite the above series of equations in matrix form A = FX where F is a 4x3 matrix containing 
the fractional elemental compositions; A and X are column matrices containing the total amount of each 
element and the quantities of each ore respectively Secondly, calculate the percent weight for weight 
(%w/w) of each element in a 10 kg mixture of the ores containing 4 kg of Ore 1, 3 kg of Ore 2 and 3 
kg of Ore 3. 


Tump to Solution 1 (see page 89) 


5.2.2 Question 2: Chemical Connectivity 

Matrices can be used to represent molecular connectivity to show, using Maths, those atoms which are 
bonded to one another. For example, ethene, or ethylene, with the atoms numbered for easy reference, 
is shown in Fig. 5.3. 



H 2 H 4 

Figure 5.3: ethene connectivity. 
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A connectivity table is shown below, Fig. 5.4, as follows: on the diagonal elements of the matrix, place 
an a (Greek alpha); on the off-diagonal elements where an atom is bonded to the indicated atom place 
a P (Greek beta); and on the off-diagonal elements where an atom is not bonded to the indicated atom 
place a 0 (zero). The connectivity table can be conveniently written in matrix form and we will use this 
type of matrix later on for carrying out some symmetry operations and also we will use some related 
matrices for molecular orbital calculations. 



C 1 

c 2 

H 1 

H 2 

H 3 

H 4 

c 1 

a 

p 

P 

P 

0 

0 

c 2 

p 

a 

0 

0 

p 

P 

H 1 

p 

0 

a 

0 

0 

0 

H 2 

p 

0 

0 

a 

0 

0 

H 3 

0 

p 

0 

0 

a 

0 

H 4 

0 

p 

0 

0 

0 

a 


a P P /3 0 0 
P a 0 0 P P 
P 0 a 0 0 0 
P 0 0 a 0 0 
0 P 0 0 a 0 
0 P 0 0 0 a 


Figure 5.4: ethene connectivity table and its corresponding matrix. 

Firstly, draw up similar connectivity matrices for the following molecules: sulfurhexafluoride (SF ), 
ethane, and diborane (B 2 H 6 ), all shown in Fig. 5.5. 


F 1 

H 1 

F 2 | F 6 

H2 ^ 

/f\ 

1 

/ 

F 3 | F 5 
F 4 

H 3 


H 4 H 1 H 3 H 5 

/ 

C 1 —C 2 

\ H 5 

Re H 2 H 4 H 6 


B 1 B 2 


Figure 5.5: sulfurhexafluoride, ethane, and diborane. 


Secondly, draw up connectivity matrices for the following unsaturated molecules and radicals; ethene 
(Fig. 5.3) and those shown in Figs. 5.6 and 5.7 but now only consider the heavy atoms and ignore the 
connectivity of the hydrogen atoms. 


C 2 H 


hLC 1 < 


C 3 H 2 H.C 1 


C 2 H- 

// 


-C 3 H 

\ 


C 4 H„ 


C 3 H= 

/ 

HC 2 

\ 

CH 


: C 4 H 

\ 


/ 

C 6 H 0 


C 5 H 


Figure 5.6: allyl radical, butadiene, and hexatriene. 


82 


Download free eBooks at bookboon.com 





















Advanced Maths for Chemists 


Week 5: Chemistry and Matrices 


HC 


C 2 H 

/©\ 

H v ~ y t 


C 3 H 


HC 2 -C 3 H 

HC 1 -C 4 H 


C 3 H : 

/ 

HC 2 

\ 

C 1 H- 


: C 4 H 

\ 


/ 

-C 6 H 


C 5 H 


Figure 5.7: cyclopropenyl radical, cyclobutadiene, and benzene. 


Tump to Solution 2 (see page 89) 

5.2.3 Question 3: Molecular Symmetry Operations, LinearTransformations 

The positions of atoms in a molecule can be defined in terms of Cartesian coordinates x,y, z. In the analysis 
of molecular geometry there is often a need to calculate atomic coordinates and then subsequently find 
their new positions (using matrices) due to the symmetry operations of reflection, inversion, rotation 
and rotation-reflection. 


Consider the anticlockwise rotation, around the z-axis, of an atom with coordinates X(x, y) through an 
angle 6 to the new position X'(x\ /) as in Fig. 5. 8, 



Figure 5.8: rotation of an atom at position X to position X’. 


What is the appropriate 2x2 matrix R for this rotational transformation such that, 

X'= XR 


Assuming that the z-axis is pointing out from the xy-plane (the right hand rule), how would matrix R 
need to be modified for a more general anticlockwise rotation about the z-axis of a point with coordinates 
(x, y, z)? 


Tump to Solution 3 (see page 92) 
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5.2.4 Question 4: Three-Dimensional Molecular Symmetry Operations 

For transformations in three dimensions, find matrix equations that represent the following operations. 

1. Anticlockwise rotation about the v-axis. 

2. Anticlockwise rotation about the y- axis. 

3. Reflection in the xy- plane. 

4. Reflection in the yz- plane. 

5. Reflection in the zx-plane. 

6. Inversion through the origin. 

Tump to Solution 4 (see page 94) 

5.2.5 Question 5: Molecular Structure Analysis of Cyclic Molecules 

Calculate a set of Cartesian coordinates with respect to the origin for firstly, the carbon atoms only of 
benzene given that the C-C bond distance is 1.40 A and the benzene ring is a regular hexagon, Fig. 5.9. 
Secondly, the carbon atoms only of the cyclopentadiene anion given that the C-C bond distance is 1.40 A 
and the five membered ring is a regular pentagon. To tackle this question find the coordinates of C 1 and 
then apply the appropriate symmetry operations to find the coordinates of the remaining carbon atoms. 
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Tump to Solution 5 (see page 95) 


5.2.6 Question 6: Methane Molecular Structure Analysis 

Calculate a set of Cartesian coordinates to 4 significant figures with respect to the origin for the carbon 
and hydrogen atoms of a methane molecule as shown in Fig. 5.10 given that the C-H bond distance is 
1.091 A and the H-C-H tetrahedral angle is 109.47°. 


For ease of calculation note that H 1 , H 2 and H 3 form an equilateral triangle and lie in a plane parallel 
to the yz-axis. Thus find the coordinates for H 4 which gives H 1 by rotation around the z-axis, and then 
find H 2 and H 3 by +120° and -120° rotation about the x-axis of H 1 . 


yi H i 


H 4 —■< 

zS 


/ 


X 


*Y n 3 

H 2 


H 




H 2 H 3 

Figure 5.10: methane molecular structure analysis 


Tump to Solution 6 (see page 97) 


5.2.7 Question 7: Ethane Molecular Structure Analysis 

Using the coordinate data from the solution to Question 6, and so find the coordinates for the carbon 
and hydrogen atoms in ethane for the eclipsed (Fig. 5.11 on the left) and staggered (Fig. 5.11 on the 
right) conformations. The C-C bond length is 1.541 A and C-H bond distance is 1.091 A as above and 
all bond angles are tetrahedral. You may find it helpful in each case to make a molecular model and 
consider the symmetry relationships between corresponding atoms in each of the structures. 


u4 Ay h 1 


X 


H 6 


\ 

C 2 - _ 

l/z Y|_| 

□ 5 U2n 


/ 

-01 


H 5 


H 2 



Figure 5.11 : ethane eclipsed and staggered conformations. 

Tump to Solution 7 (see page 98) 
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5.2.8 Question 8: Chemistry and Simultaneous Equations 

There is a linear relationship between the Celsius and Fahrenheit temperature scales such that 

T v = mT c + C 

Given that the freezing point of pure water is 32°F and 0°C, and the boiling point of pure water 212°F 
and 100°C, draw up a pair of associated linear equations, transform into matrix form and hence obtain 
values for the m and C. 


Tump to Solution 8 (see page 99) 


5.2.9 Question 9: The Speed of Chemical Reactions 

The Arrhenius equation is an important equation in reaction kinetics 

*= 

where k is the rate constant, A is the pre-exponential factor, £ a is the activation energy, T is the temperature 
in kelvin, and R is the gas constant. This can be rearranged into a linear form by taking natural logarithms 
of both sides, 

E 1 

In k= In A -* — 

R T 


If the value of k can be measured at a series of temperatures T, then the two unknown parameters A 
and £ a can be calculated by solving the resulting simultaneous equations. 

The initial stage of the reaction between ammonia and nitrous oxide follows second order kinetics. Given 
that the rate constant at 600 K is 0.385 L mol 1 s 1 and at 716 K is 16.0 L mol 1 s' 1 , calculate the activation 
energy £ a , and the pre-exponential factor A to 3 significant figures. 


Tump to Solution 9 (see page 100) 

5.2.10 Question 10: Chemical Analysis of Well Separated Spectral Regions 

For a single absorbing molecular species in solution, the Beer-Lambert Law states that 

A - sic 
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where A is the absorbance, e is the molar absorptivity (extinction coefficient), l is the path length of the 
cell containing the absorbing solution (usually 1 cm) and c is the concentration of the solution in mol L' 1 . 
If a mixture has spectra for each of the individual component compounds of the mixture which has 
regions (shown in Fig. 5.12 as dotted vertical lines) where each component dominates the spectrum, then 
we can use straight forward matrix algebra to determine the concentrations of the mixture components. 



Wavelength A 


Figure 5.12: spectra with dominant regions. 
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For a solution containing two absorbing species X and Y the overall absorbance measured (A total ) would 
the sum of the individual absorbances at wavelengths chosen so that one of the components is the 
major contributor (the black dashed lines in Fig. 5.12). For a two component mixture we need two such 
wavelengths to give two simultaneous equations 


A student in the teaching lab measured the following absorptivities using a 1 cm path length cuvette for 
the pure components known to be in the mixture 


//nm 

£ x /L mol 1 cm 1 

e Y /L mol 1 cm 1 

285 

16441 

3872 

338 

3993 

6423 


Figure 5.13: molar extension coefficients for pure X and Y. 

The student was given an unknown mixture of the two compounds X and Y and measured in a 1 cm 
cuvette the total absorbances of 0.952 at 285 nm and 0.553 at 338 nm. Find the concentrations of X and 
Y in the mixture by rewriting the data in matrix form for the two simultaneous equations and solving 
them by the appropriate matrix manipulations. 


Tump to Solution 10 (see page 101) 

5.2.11 Question 11: Chemical Analysis of the Spectra of Three or More Components 

A student in their lab class measured the following spectroscopic data for some pure samples of three 
indicators thymol blue (TB), semithymol blue (STB) and methylthymol blue (MTB), Fig. 5.14. 



TB 

STB 

MTB 

Mixture 

Mixture 

//nm 

£ x /L mol -1 cm -1 

£ y /L mol 1 cm" 1 

£ x /L mol 1 cm -1 

A 

Concn. 

455 

4800 

11100 

18900 

0.4745 

3.5230E-06 

485 

7350 

11200 

11800 

0.3394 

6.5143E-06 

545 

36400 

13900 

4450 

0.3095 

2.0385E-05 


Figure 5.14: absorptivities for pure TB, STB and MTB and absorbances for a mixture. 


Rewrite the data in matrix form and by appropriate matrix manipulation; find the concentrations of TB, 
STB and MTB in the mixture. You should attempt the manual approach which we have used up until 
now and find the difficulties involved. I will show you how to use a spreadsheet with matrix functions 
to solve this problem in section 5.3.11. 


88 


Download free eBooks at bookboon.com 



















Advanced Maths for Chemists 


Week 5: Chemistry and Matrices 



Figure 5.15: 2-dimensional model of thymol blue, S= yellow, 0 = red, H is not shown. 

Tump to Solution 11 (see page 103) 


5.3 Week 5 Tutorial Solutions 

5.3.1 Solution 1: Chemical Ore Analysis 

Firstly, the Analytical Chemists results written in matrix form are given below, where A = FX. 



^Pb 


[0.1 

0.2 

0.3 


\ 

A = 

^Fe 

F = 

0.2 

0.3 

0.3 

X = 


^Cu 


0.6 

0.2 

0.2 





0-1 

0.3 

0.2 


w 


For the 10 kg of the mixed ores, the product of a 4x3 matrix and a 3x1 matrix is a 4x1 matrix, 


^Pb 


0.1 

0.2 

0.3 


/a) 


0.1x4 + 0.2x3 + 0.3x3 


1.9 

^Fe 


0.2 

0.3 

0.3 




0.2X4+ 0.3X3+ 0.3X3 


2.6 

^Cu 


0.6 

0.2 

0.2 


J 

a 


0.6X4 + 0.2X3 + 0.2X3 


3.6 



\1 

0.3 

0 .2j 


VI 


0.1x4 + 0.3x3 + 0.2X3^ 


H 


The 10 kg of mixture has Pb = 19%, Fe = 26%, Cu = 36% and Mn = 19%. 

Return to Question 1 (see page 80) 

5.3.2 Solution 2: Chemical Connectivity 

Firstly we are asked to draw up connectivity matrices for SF 6 , ethane and diborane. We note that although 
ethane and diborane have the similar empirical formulae their connectivities are very different. 
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F 2 

F 3 

F 4 
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P 

P 

P 

P 

F 1 

p 

a 

0 

0 

0 

0 

0 

F 2 

p 

0 

a 

0 

0 

0 

0 

F 3 

p 

0 

0 

a 

0 

0 

0 

F 4 

p 

0 

0 

0 

a 

0 

0 

F 5 

p 

0 

0 

0 

0 

a 

0 

F 6 

p 

0 

0 

0 

0 

0 

a 



a P P P P P P 
Pa 0 0 0 0 0 
P 0 a 0 0 0 0 
p 0 0 a 0 0 0 
P 0 0 0 a 0 0 
p 0 0 0 0 a 0 
p 0 0 0 0 0 a 
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H 2 


H 3 


H 4 

/ 
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\ H 5 
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H 2 
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H 4 

H 5 
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a 

p 

P 

P 

P 

0 

0 

0 
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p 

a 

0 

0 

0 
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p 

P 
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p 

0 

a 

0 

0 

0 

0 

0 
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p 

0 

0 

a 

0 

0 

0 

0 
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p 

0 

0 
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a 

0 

0 
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0 

p 

0 

0 
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a 

0 

0 
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0 

p 

0 

0 

0 

0 

a 

0 
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0 

p 

0 

0 

0 

0 

0 
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a p 
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p 0 
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0 p 
0 p 
0 p 


p p 
0 0 
a 0 
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0 0 
0 0 
0 0 
0 0 


P 0 
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B 1 

B 2 

H 1 

H 2 

H 3 

H 4 

H 5 

H 6 

B 1 

a 

0 

P 

P 

P 

P 

0 

0 

B 2 

0 

a 

0 

0 

P 

P 

p 

P 

H 1 

p 
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a 
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0 
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0 
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a 

0 

0 

0 
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H 3 

p 

p 

0 

0 

a 

0 

0 

0 

H 4 

p 

p 

0 

0 

0 

a 

0 

0 

H 5 

0 

p 

0 

0 

0 

0 

a 

0 

H 6 

0 

p 

0 

0 

0 

0 

0 

a 


aO/3/3/3/300 

0a00/3/3/3/3 

pOaOOOOO 

pOOaOOOO 

PfiOOaOOO 

PPOOOaOO 

OpOOOOaO 

0/?00000a 


Secondly, we are asked to draw up connectivity matrices for the unsaturated acyclic molecules and radicals 
ethene, allyl radical, butadiene, and hexatriene; and then their cyclic counterparts cyclopropenyl radical, 
cyclobutadiene and benzene. For all these species we will consider only the carbon atoms. 
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HC 2 -C 3 H 

HC 1 -C 4 H 
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p 
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p 

0 

0 
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0 

p 

a 

p 

0 

0 

c 4 

0 
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p 
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0 
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p 
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0 

0 
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a p 0 0 Op 
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0 P a P 0 0 
0 0 P a P 0 
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P 0 0 0 p « 


Return to Question 2 (see page 81) 

5.3.3 Solution 3: Molecular Symmetry Operations, Linear Transformations 



If the distance from the origin is r then for X and X' we have, 

X (x = rcoscr y - r sin cr) X ' (pc' = rcos (a +Q) y' = rsin (a +0)) 

Expanding the double angles using the trig formula cos(A + B) = (cosA cos B) - (sinA sin B) and the 
formula for sin(A + B) = (sinA cos B) + (cosA sin B). For X' we have, 

x'= r(cosa cos0 - sin a sin0) y ,= r(sina cos0 + cosasin©) 


Noting that the change from (x, y) to (x', /) only involves the value of Q (as r and a are constants) we 
can put the change due to the rotation, i.e. the transformation, into a matrix form, 


x ? 


x cos 6 — ysin# 
^xsin0 + ycos d 



cos 6 — sin 6 W x 

^sin0 cos0 ) 
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The rotation matrix around the z-axis R is 


R = 


cos # 
sin# 


—sin# 

cos# 


Extending the rotation matrix R to the more general anticlockwise rotation about the z-axis of a point 
with coordinates (x, y, z) 


M 


xcosd 

—y sin# 

zX 0 


lx') 


COS d 

—sin# 

o' 

/ \ 
X 

y' 

— 

xsinO 

y cos 6 

zXO 


y' 

— 

sin d 

COS0 

0 

y 

\ z 'i 


{ xXO 

yXO 

zX 1 } 


\ z 'i 


0 

0 

!/ 

W 


R = 


cos# 

sin# 




0 


—sin# 

cos# 

0 


0 

0 

1 


\ 

/ 


Return to Question 3 (see page 83) 
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5.3.4 Solution 4: Three-Dimensional Molecular Symmetry Operations 



Figure 5.17: three-dimensional Cartesian coordinates. 

We have already looked at the solution for the rotation around the z axis in section 5.3.3 above, Now 
anticlockwise rotation about the v axis by an angle 9 is. 


X ' 1 


1 

0 

0 

/ \ 
X 

y’ 

= 

0 

cos d 

—sin# 

y 

U'J 


0 

sin# 

COS# j 

u 


and anticlockwise rotation about the y-axis by an angle 9 is, 


\ 

JC'I 


cos# 

0 

sin# 

/ 

X 

y ' 

— 

0 

1 

0 

y 

z 'i 


sin# 

0 

cos#^ 

\ z l 


Reflection in the xy plane only changes the z coordinate to -z 
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l 

0 

0 ' 

/ \ 
X 

y ' 
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\ z 'i 


1° 

0 

- 1 / 

\ z l 


and reflection in the yz plane only changes the x coordinate to -x 
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X ' 1 


1-1 

0 

o' 

/ \ 
X 
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\° 
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!/ 

\ z l 


and the final reflection, in the zx plane only changes the y coordinate to -y 


X') 


i 

0 

\ 

0 

1 

X 

y' 

— 

0 

-1 

0 

y 

z 'l 


1° 

0 

l l 

\ z i 


Inversion through the origin will change all three coordinates to their negative values 


lx') 


-l 

0 

0 ' 

/ 

X 

y' 

= 

0 

-1 

0 

y 

\ z 'j 


1° 

0 

- 1 / 

\ z l 


Return to Question 4 (see page 84) 
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5.3.5 Solution 5: Molecular Structure Analysis of Cyclic Molecules 

The C-C distance in benzene is 1.40 A with the angle between neighbouring C-atoms of 360°/6 = 60°. 
So neighbouring C-atoms form an isosceles triangle with the origin, so half the angle and half the bond 
distance gives us a right angle triangle with the hypotenuse of (1.40 A/2)/sin(60°/2) = 1.40 A. In Fig. 
5.18 the coordinates of C 1 are (1.40, 0, 0) after dropping the units for clarity. 



Figure 5.18: benzene carbon atoms. 


C 2 is symmetry related to C 1 as a rotation around the z-axis of 300° 


l~ 2 

C x 


cos 300 —sin 300 0 

^1.40^ 


, (l.40Xcos300)- (0Xsin300)+ (0x0)l 
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0 / 


C 3 is a reflection in the yz plane of C 2 and only changes the * coordinate of C 2 to -x 
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1 

O 
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C 4 is a reflection in the yz plane of C 1 and only changes the * coordinate of C 1 to -x 
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C 5 is an inversion through the origin of C 2 will change all three coordinates to their negative values 
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C 6 is a reflection in the zx plane of C 2 and only changes the y coordinate of C 2 to -y 


C x 


1 

0 

o' 

/ \ 

0.70 


\ 

0.70 

cV 

= 

0 

-1 

0 

-1.21 

~ 

1.21 

C 6 z 


1° 

0 

V 

\ 0 1 


0 / 
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The C-C distance in the cyclopentadiene anion is 1.40 A with the angle between neighbouring C-atoms of 
36075 = 72°. So neighbouring C-atoms form an isosceles triangle with the origin, so half the angle and half 
the bond distance gives us a right angle triangle with the hypotenuse of (1.40 A/2)/sin(72o/2) = 1.19 A. 
So in Fig. 5.19 C 1 has the coordinates (1.19, 0, 0) after dropping the units for clarity. 



Figure 5.19: cyclopentadiene anion carbon atoms. 


C 2 is obtained by an anticlockwise rotation around the z-axis of 360° - 72° = 288°, 


C 2 xl 


cos288 

—sin 288 

\ 

0 

/ \ 

1.191 


(l.l9Xcos288) — (0Xsin299) + (0x0) 


\ 

0.37 

C 2 j 

~ 

sin 288 

cos 288 

0 

0 

= 

(1.19Xsin 288) + (OXcos 288) + (0X0) 

= 

-1.13 

\c 2 zj 


l 0 

0 

V 

0 i 


(1.19X0) + (0X0) + (Ox l) J 


0 / 


C 3 is an anticlockwise rotation of C 2 by 288° around the z-axis, 


C x 


cos 288 —sin 288 0 

/ \ 

0.37 


(0.37Xcos288) + (1.13Xsin288) + (OxO)] 


\ 

[-0.96 

c 3 j 

= 

sin288 cos288 0 

-1.13 

= 

(0.37Xsin288) - (l.l3Xcos288) + (0x0) 

= 

-0.70 

\ c 7 


l 0 0 V 

\ 0 i 


(0.37X0)-(1.13X0)+ (0xl) j 


0 / 
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C 4 is a reflection in the zx plane of C 3 which only changes the y coordinate, 


C 4 x' 


1 

0 

o' 

-0.96' 


1- 0.96' 

cV 

= 

0 

-1 

0 

-0.70 

— 

0.70 

c 4 zj 


1° 

0 

V 

0 1 


0 / 


C 5 is a reflection in the zx plane of C 2 and only changes the y coordinate, 


C 5 x | 



0 

o' 

1 0.37 1 


\ 

0.37 

C 5 7 

~ 

0 

-1 

0 

-1.13 

= 

1.13 

\c 5 zj 


1° 

0 

A 

\ 0 i 


0 / 


The take home message is that without the use of matrix algebra of symmetry operations these structural 
problems in Chemistry would involve some complicated non-trivial trig calculations, I know which I 
would rather do! 


Return to Question 5 (see page 85) 


5.3.6 Solution 6: Methane Molecular Structure Analysis 

/in 1 

A 


H 4 - 

z' 


x 


H 3 


H 2 


z v\ 

H 2 H 3 


Figure 5.20: methane structure. 


In Fig. 5.20 the coordinates of C 1 are (0, 0, 0) and the coordinates of H 4 are (-1.091, 0, 0). 


H 1 can be found by a clockwise (or negative) rotation of 109.47° about the z-axis of H 4 . 


H'xl 


cos(-109.47) 

-sin (-109.47) 

\ 

0 

—1.091^ 


(-1.091 Xcos(-109.47)1 


0.3636' 

H l y 

= 

sin (-109.47) 

cos(-109.47) 

0 

0 

= 

-1.091 Xsin(-109.47) 

= 

1.0286 

\H‘zj 


0 

0 

1 

\ 0 i 


0 ) 


0 / 


H 2 and H 3 can be found by a rotation of +120° and -120°, respectively, about the v-axis of atom H 1 . 


H 2 x 



0 

0 ' 

I \ 

0.36361 


0.3636 1 


\ 

0.3636 

H 2 y 

= 

0 

cos 120 

—sin 120 

1.0286 

~ 

1.0286xcos 120 

= 

-0.5143 

H 2 z 



sin 120 

cos 120 1 

\ 0 i 


i 1.0286 X sin 120^ 


0.8908 / 


H 3 x 


\ 

10 0 1 


0.36361 


\ 

0.3636 | 


0.3636 

H 3 J 

= 

0 cos (— 120) —sin(—120) 


1.0286 

= 

1.0286Xcos(—120) 

= 

-0.5143 

\H 3 zj 


^0 sin(—120) cos(— 120) ^ 


0 i 


1 1.0286 X sin (—120)j 


[-0.8908/ 


Return to Question 6 (see page 85) 
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5.3.7 Solution 7: Ethane Molecular Structure Analysis 

H 6 




y h< 


H 5 


. 0 . 1 . 


H 6 , 


C 2 ^ 

M v 


X 




y w 


x 


H 5 


H 2 


H 3 


. 0 . 1 ... 

//Z \ I—13 

H 4 H 2hl 


Figure 5.21: ethane structure. 


The C-C bond length is 1.541 A and C-H bond length is 1.091 A and all bond angles are tetrahedral 
109.47°. For the eclipsed conformation (on the left of Fig. 5.21) the coordinates of C 1 , H 1 , H 2 , and H 3 
are similar to methane (Solution 6) except the x coordinates are increased by 1.541 A/2 = 0.7705 A. The 
coordinates of C 2 , H 4 , H 5 , and H 6 are related to C 1 , H 1 , H 2 , and FI 3 by reflection in the mirror in the yz 
plane i.e. the x coordinates change sign. 



X 

y 

z 


X 

y 

z 

c 1 

0.7705 

0.0000 

0.0000 

c 2 

-0.7705 

0.0000 

0.0000 

H 1 

1.1341 

1.0286 

0.0000 

H 4 

-1.1341 

1.0286 

0.0000 

H 2 

1.1341 

-0.5143 

0.8908 

H 5 

-1.1341 

-0.5143 

0.8908 

H 3 

1.1341 

-0.5143 

-0.8908 

H 6 

-1.1341 

-0.5143 

-0.8908 


Figure 5.22: ethane eclipsed coordinates. 


The coordinates of C 1 , C 2 , H 1 , FI 2 , and H 3 of the staggered conformation of ethane (on the right in Fig. 
5.21) are identical to those above in Fig. 5.22. The coordinates of H 4 , H 5 , and H 6 are related to H 1 , H 2 , 
and H 3 by inversion through the centre of symmetry, i.e. the coordinates of H 4 , H 5 , and H 6 are the same 
as H 1 , H 2 , and H 3 but with an inversion of their sign. 



X 

y 

z 


X 

y 

Z 

C 1 

0.7705 

0.0000 

0.0000 

c 2 

-0.7705 

0.0000 

0.0000 

H 1 

1.1341 

1.0286 

0.0000 

H 4 

-1.1341 

-1.0286 

0.0000 

H 2 

1.1341 

-0.5143 

0.8908 

H 5 

-1.1341 

0.5143 

-0.8908 

H 3 

1.1341 

-0.5143 

-0.8908 

H 6 

-1.1341 

0.5143 

0.8908 


Figure 5.23: ethane staggered coordinates. 


Return to Question 7 (see page 85) 
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5.3.8 Solution 8: Chemistry and Simultaneous Equations 

There is a linear relationship between the Celsius and Fahrenheit temperature scales that for two different 
temperatures which may then be written in matrix form and substituting for the freezing point and 
normal boiling points of water gives, 


T\ = mT l c + C 
T\= mT 2 c + C 







0 1 ^ lm\ 

Woo ljWi 


Writing the matrix equation in the standard way for simultaneous equations gives, 

0 = ( 32 

100 1/\C/ \ 212 


This is in the form A x = b and using the matrix algebra described previously, the determinant |A|, the 
matrix of cofactors C, the transpose of the matrix of cofactors C, and the inverse A -1 are, respectively, 


A 


0 

100 


= 0X1 - 1X100 =-100 
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' + 1 -iooj 

1-1 + o I 


c T = 


+1 -1 
-100 +0 



1 

-100 


+1 

l-100 


-1 
+ 0 ; 


The solution for the unknowns m and C is obtained from x = A 1 b 


/ 

l 


—!—( 1 -> 
-100 \-100 0 



1 

-100 


(lx 32) — (lx212) 
J—100X32)+ ( 0 x 212 ), 


1 

-100 


' -180 \ 
-3200/ 


1.8 

32 


T f = 1.8 T c + 32 


Return to Question 8 (see page 86) 


5.3.9 Solution 9: The Speed of Chemical Reactions 

/ln*,Wl -l/r^/ln^ 
(lnfc 2 j (l -1/T 2 ]\eJR] 


In the standard matrix form A x = b this is rearranged to 


[l -l/r i \/ln^\ 

(i -\it 2 )\ejr) 


I In k |' 
\ln^ 2 , 


Substituting in the experimental values without units for clarity, and for the moment working to extra 
significant figures to avoid rounding errors, 

'l —1/600\ / ln^ \ _ /in0.385) 

^1 -1/716/ \£ a //?j \ In 16.0/ 

(l —1.6667X 10“ 3 W ln^4 \ /—0.9545^ 

\1 -1.3966xl0' 3 /^ a //?j \ 2.7726 ) 


A | = (— 1.3966X10 3 ) — (— 1.6667X 10 3 ) = 2.701xl0~ 4 

c J+(-1.3966xloi — 1 j _ /— 1.3966X10 -3 -l] 
— (—1.6667 X10 3 ) +1/ ( + 1.6667X 10~ 3 +1/ 


-1.3966X1 O' 3 +1.6667X10 3 

-1 +1 
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Solving for the unknowns matrix x using, 


x — A b x — y—\C T b 
\A\ 


1 


1 


x — 


2.701X10 


2.701X10 

l 


-1.3966X10' 3 + 1.6667x10^ / 

-1 +1 


-0.9545 

2.7726 


1.3331X10 -3 + 4.6211X10 -3 ] 
0.9545 + 2.7726 


1 


2.701 Xl0~ 


5.9542X10 

3.7271 


-3 ' 


_ [22.044] 
* \ 13799 J 


Thus In A = 22.044 and the pre-exponential factor to three significant figures is, 

A = 3.75x 10 9 L mol' 1 s" 1 


£ a AR = 13,799 so £ a = 13,799x8.314 J mol 1 K 1 = 114.72xl0 3 J mol 1 K 1 or in the conventional units 
and to three significant figures, 

E a = 115 kJ mol -1 


Return to Question 9 (see page 86) 


5.3.10 Solution 10: Chemical Analysis of Well Separated Spectral Regions 

For the binary mixture then writing this in matrix form, entering the student’s lab results (without units 
for clarity), rearranging the matrix equation into the standard simultaneous equations form gives A x = 
b and solving this by matrix algebra 


^285 ^X285 ^ P^-] + ^ Y285 ^ 

^338 “ ^X338 ^ P^l + ^ Y338 ^ 



I 8 X2S5 
\ ^ X338 



/o.952 Wl6441 3872 Wc x ' 
\0.553 J \ 3993 6423 J ^C Y 


16441 3872 \(C X \ (o.952] 

^ 3993 6423 J \C Y j \0.553 j 
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\A\= 16441 (6423)-3872(3993)= 1.056x 10 8 - 1.546x 10 7 = 9.014x 10 7 


'+6423 -3993 \ 

-3872 +16441/ 


'+6423 —3872 \ 

{ —3993 +16441/ 


1 

II 

X 

x= —— C T /> 


\A\ 


1 / + 6423 —3872 \ /0.952\ _ 1 [ (0.952x6423) + (0.553 X—3872) \ 

9.014X 10 7 \—3993 +16441 / \0.553 ) 9.014X 10 7 \(0.952x-3993) + (0.553X 1644l )) 

1 ( 6115-2141 1 /3974| 

X 9.014X10 7 \ — 3801 + 9092 / 9.014X10 7 \5291 ) 

_l 4.4087 XlO^ 5 \ 

X _ \5.8698X10“ 5 ] 


The unknown concentrations from the well-separated spectral regions are to three significant figures 


Cx = 4.41 x 10 ~ 5 molL ~ 1 and c Y = 5.87x 10 _5 molL 
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Where there is a substantial overlap of the spectra of the pure components (Fig. 5.24) we cannot used 
matrix algebra methods to solve the simultaneous equations. They are solved by a linear least squares 
spreadsheet using the Solver algorithm (Parker 2013, section 1.1.2), (Blanco et al 1989), and (Harris 
2007, p 402). 



Figure 5.24: spectra of pure components with substantial overlap 


Return to Question 10 (see page 86) 

5.3.11 Solution 11 : Chemical Analysis of the Spectra of Three or More Components 

This is the same fundamental problem as above, extended to a ternary mixture. It is assumed that the 
concentrations of the dyes are low and their absorbances at a range of wavelengths are additive, i.e. they 
are each dominant in at a given wavelength and they don’t interact chemically with one another. First 
let us try the manual matrix approach. 

^455 “ ^TB 455^ C TB + ^STB455 ^ C STB + £ MTB 455 MTB 
^485 ” ^TB485^ C TB + ^STB485 ^ C STB + ^ MTB 485 MTB 
^545 “ ^TB545^ C TB + ^STB545 ^ C STB + ^MTB545 ^ C MTB 


Writing this in matrix form, entering the students lab results (without units for clarity) and then 
rearranging the matrix equation into the standard simultaneous equations form we have 


A 

^455 


^ TB455 ^ STB 455 ^ MTB 455 


C TB 

oo 

— 

^TB485 ^ STB485 ^MTB485 


C STB 

\ 4 54 5 


TTB545 ^STB545 MTB 545 j 


i C MTB 


\ 

0.4745 


4800 11100 18900^ 


C TB 

0.3394 

= 

7350 11200 11800 


C STB 

\0.3095 j 


^36400 13900 4450 j 


C MTB 


\ 


\ 


\ 

4800 11100 18900] 


C TB 


0.4745 

7350 11200 11800 


C STB 

= 

0.3394 

\36400 13900 4450 j 


i ^ MTB , 


^0.3095 / 
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but trying to solve this manually gives us problems , here is the determinant 


\A\ = +4800(11200x4450- 11800x 13900)- 11100(7350x 4450- 11800x 36400) 

+ 18900(7350x 13900- 11200x 36400)- 7350(11100x 4450- 18900x 13900) 
+ 11200(4800x4450- 18900x 36400)- 11800(4800x 13900- 11100x 36400) 
+36400 (lllOOx 11800- 18900x 11200)- 13900(4800x4450- 18900x 36400) 
+ 4450(4800 x 11200- lllOOx 7350) 


The above determinant is awful! The manual method of matrix algebra for this type of numerical matrix 
equation has been useful for initial learning purposes also for small 2x2 matrices and matrices with 
“sparse” numbers. The manual method is also useful for examination questions where computers are not 
allowed! However, the manual method is a problem when finding the determinant and the matrix of 
cofactors for 3x3 or larger matrices and for those with numerically “rich” matrices as in the one above. 
It is slow, clumsy and prone to errors. Lets see how easy it is to use a spreadsheet and the inversion of 
matrix A. and then directly use x = A~ l b which is shown in Fig. 5.25. 

x= A cT b X= A' l b 
\A\ 



A 

B 

C 

D 

E 

F 

G 

1 


TB 

STB 

MTB 

Mixture 

Mixture 


2 

2 /nm 

s x l L mol 1 cm" 1 

£ y /L mol 1 cm" 1 

s x /L mol 1 cm' 1 

A 

Concn. 


3 

455 

4800 

11100 

18900 

0.4745 

3.5230E-06 

<— c 

TB 

4 

485 

7350 

11200 

11800 

0.3394 

6.5143E-06 

* C STB 

5 

545 

36400 

13900 

4450 

0.3095 

2.0385E-05 

< C MTB 


Figure 5.25: matrix algebra with a spreadsheet. 


The wavelengths are not used in the calculation, they are just helpful reminders. The matrix A is in the 
array B3:D5 and the matrix b is in array E3:E5. We select the three empty cell F3:F5 and type in the 
formula window the following array formula 

=MMULT (MIN VERSE(B 3 :D 5),E3 :E5) 


Do not press enter or return as it is an array function rather than a simple Maths function. Instead on 
a PC press Shiff+Control+Return or on Mac Command+Return and the solution appear in the array 
F3:F5. The spreadsheet commands are MMULT for matrix multiplication the two matrices. The first 
matrix is the inverse of A (i.e. A -1 ) MINVERSE and the second is the absorbances of the mixture, b. 


c ra = 3.52x 10~ 6 mol L~ 1 c STB = 6.51x 10 _6 molL _1 c MTB = 2.04x 10“ 5 molL _1 


Return to Question 11 (see page 88 ) 
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6 Week 6: Chemistry, Matrices and 
Huckel Theory 

Understanding the spectroscopy, electronic properties and the chemical reactions of conjugated n- systems 
is important for biological, medicinal and industrial reasons. A “simple” approach is Huckel molecular 
orbital theory where the molecular orbital energies and wavefunctions are obtained by using matrix 
algebra. I am using the example of ethene to explain in some detail the process, as often there are some 
steps omitted when presented in lectures or textbooks which don’t help our understanding of what is 
happening. Hang in there it is much easier to solve them in practice but I am starting by giving the 
background quantum mechanics. 
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Finding the energies (eigenvalues) and the wavefunctions (eigenfunctions) of the allowed quantum states 
turns up frequently in Chemistry. One of the first approaches to investigating the electronic states of 
unsaturated organic molecules with conjugated 7 r-systems is the Hiickel MO (HMO) method. The Hiickel 
approximation only considers the 7T-orbitals and ignores the a-orbitals, the o skeleton is already assumed 
in the HMO approach. The reason for spending so much time on this “qualitative” theory as opposed 
to the more sophisticated and computer based methods of ab initio or semi-empirical theories is that 
we can really get at the “heart” of Hiickel theory, otherwise quantum chemistry will be much more 
challenging and more of a “black box” for most people. Hiickel theory also gives us a clear insight into 
the molecular orbitals of a very important range of molecules. 

Using ethene as example I will work through it to show the details. The wavefunction ¥ for the 7 T-orbitals 
are made up by a linear combination of the two 2 p atomic orbitals (AOs) <j) l and 0 2 of the two carbon 
atoms contributing to the 7T-bond. 


W — c l (j) l + c 2 (p 2 

where c ] and c 2 tell us the relative contributions of these two p AOs. Different wavefunctions can have 
different values for these coefficients. The Hiickel method starts with the Schrodinger equation, 

HW = EW 

where H is the Hamiltonian or energy operator applied to the wavefunction ¥ and E the energy of the 
molecular orbital. As we shall soon see, Hiickel method expresses the energy E as the sum of two terms 
called a (the energy of an electron in a 2 p-orbital) and /? (the interaction energy between two neighbouring 
p orbitals). Both a and /? are negative quantities and at the moment unknown but importantly have 
become independent of any particular molecule. Therefore we can compare different molecules with 
one another. Substituting for ¥ gives 

c x H (p x + c 2 H (j) 2 — c l E(j) l + c 2 E (p 2 

If this is multiplied by <f)* (the complex conjugate of (f)j) and integrated over all space, and separately if 
the equation above is multiplied by <j 6 2 * and integrated we get a pair of simultaneous equation for ethene, 

c 1 J H(/) l dv + c 2 J (j) x * H (j) 2 dv = c x E J (j) x * (j) x &v + c 2 E § (j)*(j) 2 dv 
c x J 0 2 *//0 1 dv+ cj (j) 2 * H (j) 2 dv — c x E J (p 2 * (j) x dv + c 2 E J 0 2 *0 2 dv 

If the Hamiltonian integrals (on the left hand side) are written H.. for shorthand and the overlap integral 
(on the right hand side) are written as S these two simultaneous equations may be rewritten as, 

c x H xx + c 2 H 12 = c x E S xx + c 2 E S X2 
c x H 2X "I - c 2 H 22 c x E S 2X "I - c 2 E S 22 
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This pair of simultaneous equations may be rearranged to a clearer form 

(H n ~ES n )c l + (H n -ES l2 )c 2 = 0 
(H 2 1 ~ ES 21 )c x + (H 12 ~~ ES 22 )c 2 = 0 

Huckel then made the following approximations in order to be able solve the simultaneous equations 
by hand 


H 22 G 

*^11 “ $22 ~ 1 


H 12 = H 2 \ ~ P 

^ 12 “ *$21 “ 0 


and H = 0 for any non-neighbouring 7r-carbon atoms (which are not present in ethene) and then writing 
the pair of simultaneous equations in matrix form gives us the Huckel matrix equation. 


a 

i 

H=(° 

The Huckel Matrix 

( p a-EJ 

W loj 



The first matrix on the left looks very similar to our connectivity matrices developed earlier (sections 
5.2.2 and 5.3.2) but with a in the connectivity matrix replaced by a - E in the Huckel matrix. For non¬ 
trivial solutions of the Huckel matrix, the corresponding Huckel secular determinant must be equal to 0. 


a—E 

ft 

ft 

a — E 


The Huckel Secular Determinant 


To find the eigenvalues (allowed energies) we must solve the Huckel secular determinant for ethene. It 
is mathematically convenient to divide each element of the Huckel determinant by /$ 


a—E 

~T~ 

i 


i 


fi 


Making the substitution x = (a - E)//3 gives us a very simple determinant to solve for the ethene 
eigenvalues, 


x 

1 


1 

x 


= 0 


x 2 —1=0 x = ±1 


a — E 

~T~ 


= ±i 


E x — a + /3 and E 2 = a — ft The HOMO-LUMO energy gap is ^E x — E 2 ) = 2/3 
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Remembering that both a and fi are negative values, the highest occupied molecular orbital (HOMO) has 
an energy of E 1 = a + /3 for v = -1 and is n bonding. The lowest unoccupied molecular orbital (LUMO) 
has an eigenvalue E 2 = a - for v = +1 and is n * anti-bonding. The MOs are labelled with a subscript 
1 , 2, ... with increasing energy. 


To find the eigenvectors (wavefunctions) we must solve the Hiickel matrix for ethene for each of the 
allowed energies (eigenvalues). Substituting for the energy of the HOMO E = a + /3 gives (a - E ) = -/?, 


a-E x fi \ Cij _ /o\ 

l P a - E lj\ C 2l \ 0 / 


(-P 0 W^Wol ~I3c x +I3c 2 = 0 

\P _ /?/\c 2 j \0/ +/3c x - /3 c 2 = 0 


The pair of simultaneous equations implies that c x - c r The coefficients are normalized (that is the total 
probability of finding the 7T-electron in the HOMO orbital must be unity) by making the sum of the 
coefficients squared equal to unity c 2 + c 2 + ... c n 2 = 1. 


2 

C l 


+ C 1 


— 1 and c l = c 2 


thus C x — C 2 



The HOMO eigenvector (wavefunction) for ethene and its eigenvalue (energy) from Hiickel theory is, 


= (^1 + 02 ) 


E x — a + /? 
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Substituting the LUMO eigenvalue E 2 = a - (3 then (a - E ) = j3 into the Huckel matrix gives 


a — E 2 /? \i c i 


13 a — E 


2 I \ C 2] 


2 


p p 
p p 


\ C 2] 


0 j /3 c x + (3 c 2 — 0 

0 J f} c l + (3 c 2 — 0 


which means that c x - -c 2 and normalizing the coefficients 


c \ + c 2 — 1 and c x — — c 2 thus c, = -y= and c 2 = — 


The LUMO eigenvector (wavefunction) for ethene and its eigenvalue (energy) from Hiickel theory is, 

V 2 = -j^(<Pi-<P 2 ) E 2 = a — P 

The MO diagrams for the Hiickel approximation for ethene are shown in Fig. 6.1. 

| c | c 

Figure 6.1 : HOMO and LUMO of ethene. 


where the red volume are the positive-phase lobes and the blue the negative-phase lobes of the n bonds. 
On the left is the n bond and on the right the n anti-bond. The MO diagrams will use the same colour 
convention in this book. The MO energy diagram is shown below, 

E- a . 

Figure 6.2: Huckel MO energy diagram for ethene. 

For the 2 7r-electrons of ethene their total energy is 2{a + /?) = (2 a + 2/3) which compared with 2 isolated 
2p-electrons of energy 2a, means that Huckel theory shows that ethene is stabilized by 2/L 
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6.1 Week 6 Tutorial Questions 

6.1.1 Question 1: HMO Theory for the Allyl Cation, Radical, and Anion 


' C 2ne-'‘ C ' C in/ 1 ” 


Figure 6.3: allyl cation, radical and anion. 


(a) Write down the Hiickel secular determinant for the acyclic three C-atoms in terms of v, where v = 
(i a - E)/ ft. (b) Solve the secular determinant for * and then rewrite the solution in terms of (a - E)/ ft 
and so find the three allowed quantized energies (eigenvalues), (c) Evaluate the normalized wavefunction 
(eigenvector) for each energy (eigenvalue), (d) Which of the wavefunctions is the HOMO, LUMO or 
SOMO (singly occupied MO) as appropriate, for each of the allyl cation, radical and anion molecules? 


Tump to Solution 1 (see page 112) 


6.1.2 Question 2: HMO Theory for Cyclopropenyl Cation, Radical and Anion 


2 C 1c 2 C IQ 2 Q 

^z ^ ^z 

3 C 3 C 3 C 

2 tt e~ 3 tt e - 4ue' 

Figure 6.4: cyclopropenyl cation, radical and anion. 


(a) Write down the Hiickel secular determinant for the cyclic three C-atoms in terms of v, where x= (a 
- E)/p. (b) Solve the secular determinant for v and then rewrite the solution in terms of (a - E)/j3 and 
so find the three allowed quantized energies (eigenvalues), (c) Evaluate the normalized wavefunction 
(eigenvector) for each energy (eigenvalues), (d) Which of the wavefunctions is the HOMO, LUMO or 
SOMO as appropriate, for each of the cyclopropenyl cation, radical and anion molecules? 

Tump to Solution 2 (see page 115) 


6.1.3 Question 3: HMO Theory for 1,3-Butadiene 


JC 1C 

'c b 


Figure 6.5: butadiene. 


110 


Download free eBooks at bookboon.com 







Advanced Maths for Chemists 


Week 6: Chemistry, Matrices and Hikkel Theory 


(a) Write down the Hiickel secular determinant for the four C-atoms of butadiene in terms of where % 
= (a - E)//3. (b) Solve the secular determinant for v and then rewrite the solutions in terms of (a - E)/(3 
and so find the four allowed quantized energies (eigenvalues), (c) Evaluate the normalized wavefunction 
(eigenvector) for each energy (eigenvalues), (d) Which of the wavefunctions is the HOMO or LUMO? 

Tump to Solution 3 (see page 118) 


6.1.4 Question 4: HMO Theory for Cyclobutadiene 

1 C^= 2 C 

4 C^= 3 C 

Figure 6.6: cyclobutadiene. 

(a) Write down the Hiickel secular determinant for the four C-atoms of cyclobutadiene in terms of v, 
where x = (a - E)/ ft. (b) Solve the secular determinant for v and then rewrite the solution in terms 
of (a - E)/ ft and so find the allowed quantized energies (eigenvalues), (c) Evaluate the normalized 
wavefunction (eigenvector) for each energy (eigenvalues), (d) Which of the wavefunctions is the HOMO, 
a non-bonding NBMO or a LUMO for cyclobutadiene? 


Tump to Solution 4 (see page 122) 
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6.2 Week 6Tutorial Solutions 

6.2.1 Solution 1: HMO Theory for the Allyl Cation, Radical, and Anion 


x 1 0 

lxl 
0 1 x 


= 0 


The solution of the secular determinant is, 


c(x 1 2 - 1)— l(x—0)+0(l—0)= 0 x 3 — 2 x— 0 


x(x 2 — 2 ) = 0 x = 0, ± V2 


There are 3 energy levels which have the eigenvalues below, 

x = 0 x = +V2 x=—V2 

E — a E— [a — ^2 ft) E = {a+^2/3 


Since a and /3 are both have negative values so with increasing energy the eigenvalues are the following 


= a + V 2 P E 2 — a E 2 — a — V2 /? 


To find the three eigenvectors (wavefunctions) corresponding to these three eigenvalues we must 
separately solve the Hiickel matrix for allyl for each energy. For E 1 = (a + V 2 /?) then (a - E x ) = - V 2 /? 


0 

1 

8 


\ 

Cl 


(o) 


0 

Qd 

Qd 

1 


C \ 


fo\ 

13 a — E x P 


c 2 

= 

0 


P -42 p p 


C 2 

- 

0 

1 

8 

Qd 

O 


N 


H 


1 

Qd 

O 




\°l 


Which gives us the three simultaneous equations below when the matrices are multiplied out, 

— V2 P c x + (3 c 2 — 0 
P c x - 4 2 P c 2 + P c 3 = 0 
P c 2 — 'J 2 p c 3 = 0 

Subtracting the third equation from the first gives c x - c . From the first equation c 2 = V 2 c x and from 
the third equation c 2 - V 2 c 3 which when we normalize the coefficients using c x + c 2 + c 2 = 1 gives 

1 1 1 

2 V 2 2 


V 1 = 


L± . . ^3 

2 42 2 


E x =a + V2 fi 
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Figure 6.7: ^ of the allyl system (not to scale). 


Substituting into the Hiickel matrix for allyl with E 2 - a then (a - E 2 ) = 0 


a — E 2 

p 

0 1 


c i 

p 

a — E 2 

P 



0 

P 

a-E 2 j 


! C 3, 


\ 


\ 

0 p 0 


C 1 

pop 


c 2 

\0 P 0 j 


i C3 / 


Pc 2 — 0 
P c x + P c 3 — 0 
Pc 2 = 0 


From the first and third simultaneous equations c 2 = 0, and from second equation, c ] = -c 3 
If the coefficients are normalized, then c 2 + c 2 + c 2 = 1 then, 


c \ 



c 2 - 0 




r 3 

V2 


E 2 = a 



Figure 6.8: ip 2 of the allyl system (not to scale). 

Substituting into the Hiickel matrix for allyl with E 3 = a - V2 ft then (a - E 3 ) = V2 ft 


a — E 3 

p 

0 


C \ 

p 

a —E 3 

p 


C 2 

0 

p ' 

a-EJ 




'42/3 /3 0 1 


C 1 

P V2 p /3 


C 2 

0 P ^2 P) 


\° 3 I 


V2 /3 c 1 + /3 c 2 — 0 
/3 2 /3 c 2 + /3 c 3 = 0 

P c 2 + V2 /? c 3 = 0 
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Subtracting the third from the first simultaneous equations c 1 = c 3 ; from the first equation c 2 = - V2 c x ; from 
the third equation c 2 = -V2 c 3 thus when the coefficients are normalised, then c 2 + c 2 + c 2 = 1 we have, 


= 


2 



c 3 = 


2 


V' 3 = 


01 _ ^2 + ^3 

2 _ V 2 2 


£ 3 = a - V2 /3 



Figure 6.9: i/; 3 of the allyl system (not to scale). 
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Fig. 6.10 shows the MO energy diagrams for the three allyl species. 


c _ /o 0 cation 



Figure 6.10: MO energy diagram 


radical 

anion 

”1” 

-H- 

3k 

-H- 

for allyl cation, radical and anion. 


In the allyl cation CH 2 =CF1-CH 2 + with 2 7r-electrons, then \f/ ] is the HOMO and \f/ 2 is the LUMO. In 
the allyl radical CH 2 =CH-CH 2 « with 3 7T-electrons, then i// ] has 2 7r-electrons and \j/ 2 is the SUMO with 
1 7r-electron and f 3 is the LUMO. In the allyl anion CH 2 =CH-CH 2 “ with 4 7T-electrons, then \j/ ] has 2 
7 r-electrons and f 2 is the HOMO with 2 7r-electrons and f 3 is the LUMO. 


Return to Question 1 (see page 110) 


6.2.2 Solution 2: HMO Theory for Cyclopropenyl Cation, Radical and Anion 

The cyclic species differs from the linear three carbon species in that each ring C-atom is bonded to two 
other C-atoms. Thus by inspection, put x (a - E)//3 on the diagonal elements and 1 for the elements 
representing connected atoms in the secular determinant 


x 1 
1 x 
1 1 


1 

1 


= 0 


x 


x(x 2 - 1 ) - 1 (x l) + l(l-x)= 0 


x 3 -x-x +1 +1-x = 0 x 3 - 3 x+2 = 0 


This is a cubic equation and by inspection x = +1 is a root of the cubic, hence (x-1) is a factor of the 
cubic equation. We may find the other factor, which should be a quadratic, by long division and then 
find the roots of the quadratic using the quadratic formula. 

x 2 + x — 2 
* — 1 |x 3 — 3x + 2 

3 2 

X — X 

x 2 — 3x 
2 

X — X 

—2x + 2 
—2x — 2 

0 + 0 
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The quadratic factor ( x 2 +x-2 ) = 0 is solved as (x-l)(x+2) = 0 or if you prefer by using the quadratic 
formula. The cyclopropenyl molecules each have three eigenvalues of which two are identical (degenerate) 
and in increasing order of energy they are 


_ a-E _ 
x ~ p " 
E x = a+2(3 


2 


x = 


a-E 

p 


= +i 


E 2 = = a- 1 3 


_ a-E _ 

x ~ p ~ 
E 3 = a-(3 


+ 1 


The eigenvectors (wavefunctions) are obtained by solving the Hiickel matrix for cyclopropenyl molecules. 
For E 1 = a +2/1, then (a - E= -2/3, 


\ 

- 2 y 3 p \ 



<N 

1 


C 2 

<N 

1 


3 


—2[3 c x + P c 2 + f3 c ^ — 0 
/3 C; — 2 (3 c 2 + 13 c 3 = 0 
(3 c x + f3 c 2 — 2 13 c 3 = 0 


From these simultaneous equations 2c - c 2 + c } , 2c 2 - c + c 3 , 2c 3 = c 1 + c, and so we must have c } = c 2 
= c 3 , and normalizing the coefficients c 2 + c 2 + c } 2 = 1 gives, 



Vh = 


$1 $_ 2 _ $3 

V 3 V 3 V 3 


E x — a + 2 /3 



In Fig. 6.11 we see that the HMO f l is a delocalized 7r-bond involving all three C-atoms. The degenerate 
wavefunctions y/ 2 and i j/ 3 are obtained from solving the Hiickel matrix for E 2 = (a - (3) and also E 3 = (a 
- (3) so for (a - E 2 ) = [3 and (a - E 3 ) = (3 


0 

p 

\e 


p 

p 

p 


p\ 

p 

e 


= 0 


(3 c x + (3 c 2 + (3 c 2 — 0 

/3 c x + /3 c 2 + /3 c 3 — 0 

(3 c x + (3 c 2 + (3 c 2 — 0 
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The two degenerate MOs do not give independent solutions for the coefficients and we are free to choose 
the first solution as long as it satisfies the secular equation c x + c 2 + c 3 = 0 above. For y/ 2 we can choose 
that c 3 - 0 and then c x - -c 2 after normalizing the coefficients c 2 + c 2 + c 3 2 = 1 we have c x = -c 2 = iN2 


, $1 ^2 r n 

n~Ti El 



Figure 6.12: ip 2 of the cyclopropenyl system (not to scale). 


For i j/ 2 we chose that c x = -c 2 then the other degenerate wavefunction i j/ must have c x = +c 2 so substituting 
this into the secular equation c x + c 2 + c 3 = 0 gives c 3 = -2c x = -2 c r After normalizing the coefficients c 2 
+ c 2 + c 2 = 1 we have the y/ 3 wavefunction with c x = c 2 = 1/V6 and c 3 = -2/V6 and so, 


^ 3 =^= + 


2^3 


V6 V6 V6 


E 3 = a — 
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Figure 6.13: of the cyclopropenyl system (not to scale). 

Fig. 6.14 shows the MO energy diagrams for the three cyclopropenyl molecules. 

cation radical anion 

E 2 = E 3 = a - (3 -— | -— | - j— 

E=a . 

Ei=a+2/3 -H" -ff -ff 

Figure 6.14: MO energy diagram for cyclopropenyl cation, radical and anion (not to scale). 

The cyclopropenyl cation with two 2 p electrons with a 7T-electron energy of (2a+4/3) which compared 
with two isolated 2 p electrons of energy 2a has a stabilization energy of (2a+4/3) - 2a = 4/3. As we will 
see it is the most stable of the three molecules also its 7T-electrons are delocalized (aromatic) agreeing 
with Huckels (4n+2) rule for monocyclic compounds. 

The cyclopropenyl radical with three 2 p electrons with a 7r-electrons energy of 2(a+2/3) + (a -/3) = 
(3a +3/3) which compared with three isolated 2 p electrons of energy 3a has a stabilization energy of 
(3a +3/3) - 3a = 3/3. 

The cyclopropenyl anion with four 2 p electrons with with a 7T-electron energy of 2(a+2/3) + 2(a-/3) = 
(4a+2/3) which compared with four isolated 2 p electrons of energy 4a has a stabilization energy of only 
(4a+2/3)-4a = 2/3. It is also a diradical in order to satisfy Hunds rule of maximum multiplicity. 

Return to Question 2 (see page 110) 


6.2.3 Solution 3: HMO Theory for 1,3-Butadiene 


The solution of the Hiickel determinant for butadiene is shown below. 


x 1 
1 x 
0 1 
0 0 


0 0 
1 0 
x 1 
1 x 



X 

1 

0 


1 

1 

0 

X 

1 

X 

1 

~ 

0 

X 

1 


0 

1 

X 


0 

1 

X 


x 2 (x 2 - 1 )- x(x- 0)- (x 2 - 1 )= 0 


x 4 - 3 x 2 + 1 = 0 


118 


Download free eBooks at bookboon.com 












Advanced Maths for Chemists 


Week 6: Chemistry, Matrices and Hikkel Theory 


To solve this quartic equation let us substitute y = x 2 and then use the quadratic equation to solve for y 
(using an excess of significant figures to reduce rounding errors) we have, 


/-3y + l = 0 


y = 


3±V3 2 - 4X1X1 
2X1 


y = +2.61803 
x=+1.618 x = - 1.618 


7 = +0.38197 
x— +0.618 x— -0.618 


In increasing energy order the following eigenvalues (a and /} are both negative quantities) are obtained 
E x = a + 1.618 (3 E 2 — a +0.618 (3 E 3 = a - 0.618 /3 E 4 =a - 1.618)8 


E 4 = a- 1.618/3 - 


E 3 = ° ~ 0.618 P - LUMO 


E=a 


E 2 = a + 0.618 p 
E 1 = a+ 1.618/3 


-fH 

-th 


HOMO 


Figure 6.15: MO diagram for butadiene. 


Although solving this particular 4x4 secular determinant gives a quartic equation which may be solved 
manually to find the eigenvalues, as we have done above, for larger matrices this is too cumbersome. Also 
carrying out the matrix multiplication tp find the eigenvectors is less easy even for a 4x4 matrix. Secular 
determinants (the values of the column matrix on the right hand side being all equal to zero) cannot be 
solved using the spreadsheet method previously used in Week 5. Luckily we can generalize the results that 
we already have which gives us a method that is a distant relative of the particle in a box quantized systems. 

For an acyclic conjugated hydrocarbon containing n carbon atoms in the chain, 



Figure 6.16: acyclic conjugated hydrocarbon of n carbon atoms in the chain. 


The secular determinants all have the same form, 


X 

1 

0 . 

. 0 

0 

0 

1 

X 

1 . 

. 0 

0 

0 

0 

1 

x . 

. 0 

0 

0 

0 

0 

0 . 

. 1 

0 

0 

0 

0 

0 . 

. X 

1 

0 

0 

0 

0 . 

. 1 

X 

1 

0 

0 

0 . 

. 0 

1 
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for which the eigenvalues (allowed energies) are, 

£ = “ + 2 /)cos (_!tt) 

where n = number of carbon atoms in the conjugated chain; k = quantum number identifying the MO 
k = 1, 2, ... n with the cosine in radians. 


E x - a + 2 (3 cos - 
E x — a + 2 /3 X0.809 


Z7 ^0/2 71X2 

E 2 — a + 2 [3 cos—— 
E 2 — a + 2/3 X0.309 


jtX3 

E 3 — a + 2 [3 cos —— 
E 3 = a + 2 /?x(-0.309) 


Z7 ^0/9 ^><4 

E A — a + 2 (3 cos —— 
£ 4 = a + 2/? x(—0.809) 


E x — a + 1.618)8 E 2 — a + 0.618/? £3 = a — 0.618/? E A — a — 1.618/? 


These energies agree with our manual solutions found earlier. However, this generalized method becomes 
really useful for finding the wavefunctions (eigenvectors) for these larger Hiickel matrices. We have 
four different Hiickel matrices to solve, one for each energy. The coefficients are given by the following 
expression, note the proportionality sign as the coefficients have still to be normalized. 


c s oc sin 


jt ks 
n + 1 
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where k and n have the same meaning as before and s = atom number (1, 2, 3, ... n) with the sine in 
radians. 


V'l 

. JtX 1 x 1 


c i oc sin- 

S 

. JtXlX2 
c 2 oc sin--- 

. jtXlX3 
c 3 oc sin--- 

. JtX 1x4 
c 4 oc sm--- 


c 1 oc sin- 


ip 2 

. JtX2Xl 

. jtX2x2 
5 

. jtX2x3 
5 

. jtX2x4 


V' 3 

. JCX3X1 
c l oc sm--- 

. JCX3X2 
c 2 oc sm--- 

. JtX3x3 
c 3 oc sm--- 

. jtX3x4 


^4 

. jtX4X 1 


c 1 oc sm- ^ 

. JCX4X2 
c 9 oc sm- 

5 

. jtX4x3 
c 3 oc sm--- 

. jtX4x4 
c 4 oc sm--- 


V'l 


^2 ^3 

OC+0.9511 qoc +0.9511 
c 9 oc +0.9511 c 9 oc +0.5878 c 9 oc -0.5878 


c x oc +0.5878 


c 3 oc +0.9511 
c 4 oc +0.5878 


c 3 oc —0.5878 
c A oc -0.9511 


c 3 oc —0.5878 
c 4 oc +0.9511 


^4 

oc +0.9511 
c 2 oc -0.9511 
c 3 oc +0.9511 
c A oc -0.5878 


Normalizing the coefficients can be done in one step for as all four wavefunctions each have two 
coefficients squared of 0.58 78 2 and two of 0.95 ll 2 , 


/ 0.5878 \ 2 / 0.9511 V 
\ N / + \ N ) 


0.9511 V 
N ) 


' 0.5878 V 1 

[ n r 


N 2 = 2.5002 N = 1.5812 


Dividing each of the non-normalized coefficients byN = 1.5812 gives us the normalized coefficients below. 


c i 

c 2 

C 3 

C 4 


H> 1 ^2 

+0.3717 c, = +0.6015 c, 
+ 0.6015 c 2 =+0.3717 c 2 
+0.6015 c 3 =-0.3717 c 3 
+ 0.3717 c 4 =-0.6015 c 4 


^3 

+0.6015 Cj 
-0.3717 c 2 
0.3717 c 3 
+ 0.6015 c 4 


^4 

+0.6015 

-0.6015 

+0.6015 

-0.3717 




4 C 



Figure 6.17: from the top, i p, i/j 2 i p 3 ( j) t of butadiene looking down onto the 2p AOs (not to scale). 
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Fig. 6.17 illustrates the following, firstly, f l (HOMO-1) has the four 7T-electrons delocalized throughout 
the molecule. Secondly, f 2 (HOMO) is similar to our conventional drawing of two isolated 7T-bonds 
starting at carbon atoms 1 and 3 as in the full name 1,3-butadiene. Thirdly, \j/ 3 (LUMO) has a single 
isolated double bond between atoms 2 and 3. Finally, \j/ 4 (LUMO+1) has four isolated 2 p atomic orbitals 
on the carbon atoms and is completely 7r-antibonding. The ground electronic state has the four 2 p 
electrons occupying \f/ ] and f 2 with spins paired, see Fig. 6.15. 


Return to Question 3 (see page 110) 


6.2.4 Solution 4: HMO Theory for Cyclobutadiene 

We will apply the general method for cyclic molecules to cyclobutadiene. The secular determinant is, 


x 1 0 

1 x 1 

0 1 x 

1 0 1 


1 



x 


The determinant has an extra “1” in the (4,1) and the (1,4) corners compared to the acyclic butadiene, 
this is because of its “closed-up” ring structure. The general method for the energies (eigenvalues) of a 
cyclic conjugated hydrocarbon is slightly different from the acyclic version due to the slightly different 
secular determinant with the extra “1” in the (n, 1) and (1 ,n) corners. The cyclic molecule solutions for 
the eigenvalues and eigenvectors are distantly related to the motion of a particle on a ring (Parker 2012, 
sections 4.1.3 and 4.2.3). The particle may travel clockwise or anticlockwise around the ring, thus there 
are two possible solutions for the eigenvectors. 


x 1 0 ... 0 0 1 

1 x 1 ... 0 0 0 

0 1 x ... 0 0 0 

0 0 0 ... 1 0 0 

0 0 0 ... x 1 0 

0 0 0 ... 1 x 1 

1 0 0 ... 0 1 x 


= 0 


k E - a + 2/3 cos 


2 jt k 
i n 


The cosine is in radians where n = total number of carbon atoms forming the ring and k = quantum 
number identifying the MO where k may take the values, 
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k = 0,±1,±2,-,±|^| for n =4,6, 8 ,- 


Jfc = 0,±l,±2,-,± 


n — 1 


for ^ = 3,5,7, 


Note that for “even” rings k runs to the positive value of n/2 but for “odd” rings that k runs to plus and 
minus (n-l)/2. For cyclobutadiene n = 4 and so k may take the values k = 0, ±1, and +2 only i.e. there 
are four solution to the quartic equation in x derived from the 4x4 secular determinant. 


o i 7 , q 2 jt 0 
E — a + 2 3 cos —-— 

4 

°E =a+2/3[l) 


1 77 , ni ry 2 JT 1 

E — a + 2 3 cos —— 
4 

l E = a + 2/3[0) 


-1 


2n —1 

E = a + 26 cos- i ;—- 

4 

~'E = a +2/?(o) 


2 77 , a n 2 JT 2 

E — a + 2 3 cos —-— 

4 

2 E = a +2/5 (—1) 


Relabelling the energies in the convention way of increasing energy where a and are negative quantities, 


E x — a + 2l3 E 2 — a E 3 = a E^ — a—l/d 
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So there is one bonding, two degenerate non-bonding and one antibonding MOs for cyclobutadiene. 
There are four 7T-electrons, two of them are paired up in the f 1 bonding MO, the other two are placed 
one each in the two degenerate non-bonding MOs with their spins parallel because of Hund’s rule of 
maximum multiplicity. 


E 4 = a - 2 (3 - lij 4 antibonding LUMO 


E 2 


= £3 = “-) — f- 

= a+ 2P -ff 


<// 2 & <// 3 nonbonding SOMO 
(// 1 bonding 


Figure 6.18: Huckel MO diagram for cyclobutadiene. 


The 7T-electron energy is 2(a + 2/3) + 2 a = (4a + 2/3). The energy of two isolated double bonds (as in 
two ethene molecules) is 2(2a + 2/3) = (4a + 4/3). The net stabilization energy of cyclobutadiene is -2/3 
i.e. cyclobutadiene is destabilized by 2/3 compared with two isolated double bonds and thus is very 
reactive and will react as a di-radical. It is unstable above 35 K and dimerizes by a Diels-Alder reaction. 
As expected cyclobutadiene is not aromatic and it does not have aromatic stability as it does not have 
(4n+2) 7r-electrons as in Huckels rule for monocyclic compounds. 

To determine the eigenvectors for a cyclic molecule we need to solve the general Huckel matrix equation 
shown below for a generalized cyclic molecules. 
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The general solution for the wavefunctions are complex (contain an imaginary part due to motion of the 
particle in either direction around the ring) but the real functions maybe found from linear combinations 
the solutions in pairs. The two real solutions are given below, 


k+ lp CC V cos 

2 jt ks j ^ 

1 

8 

M 

sin 

1 2n ks V 

* ' 

1 n r s 

s 


n r s 


The coefficients will have to be normalized by dividing by the normalization constant N. For k = 0 are, 
0 2jtX0x1 , 2jtX0x2 , 2jtX0x 3 , 2jtX0x4 , 

Ip OC cos-—- (j) x + COS -—- 0 2 + COS -—- 03 +COS --- 0 4 

°lp OC +101 +10 2 + 103 + 104 
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The four coefficients are all in phase and form a delocalized 7r-bond around all four carbon atoms. This 
is the lowest energy Hiickel MO conventionally labelled \f/ ] and the normalization constant is N 2 = 4 so 
N= 2. 



Figure 6.19: ^ of cyclobutadiene looking down onto the 2 p AOs. 


The eigenvectors for k = +1 and k = -1 are, 

/ 

0 1 + COS 

\ 

/ 

0 l + sin 

\ 

+ V QC 001 1 0 2 +003 +1 04 

°c +10! +O0 2 -103 +O0 4 

Thus for k = ±1 the two real wavefunctions correspond to the two degenerate nonbonding SOMO 
wavefunctions. For degenerate wavefunctions then any linear combination gives a valid wavefunction 
(Levine 2009, p. 643). The sum and the difference of +1 f and _1 i// are f 2 = +1 f + _1 i// and f 3 = +1 f - ~ l f. 
They given below after normalization by diving by 2. These degenerate pair of wavefunctions each consists 
of two isolated double bonds, 



2 C 


3 C 


Figure 6.20: ip 2 SOMO of cyclobutadiene looking down on the 2 p AOs. 



Ijj OC cos 
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ip oc sin 


2jtXlXl 
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2jtX1 X2 


+ cos 


+ sin 


2jcX 1X3 
4 y 
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2 JtX 1X4 


\ * 
2jtXlX4^ 
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^3 = 


%_ <P2 + Vl + <P4 
2 2 2 2 


E 3 = a 



Figure 6.21 : ip 3 SOMO of cyclobutadiene looking down on the 2 p AOs. 


The eigenvector for k = 2 is 


2 . / 2:tX2Xl | , | 2 jtX2x2 | , 

7 /; oc cos I - - - I 0 1 +cosl -- -I (j) 2 + cos 


ip oc 


•1^1 


+i< 


2jtX2x3 , 2jtX2x4 \ , 

- 4 -j 03 + cos^- - -j0 4 

-103 + 10 4 


Thus for k = 2 the four coefficients are all out of phase with both of their neighbours forming a completely 
7r-antibonding MO. This is the highest energy Hiickel MO f 4 and the normalization constant is again N= 2. 



E 4 = a - 2/3 


2 C 


3 C 


Figure 6.22: ip 4 LUMO of cyclobutadiene looking down on the 2 p AOs. 

Although the Hiickel method is now replaced by computer based semi-empirical, ab initio , or DFT 
methods the Hiickel approximation is still great for allowing us to manually solve molecular orbitals 
and see what is behind these more sophisticated methods as well as giving us a good qualitative feel for 
the properties of the important class of both acyclic and cyclic conjugated molecules. 

Return to Question 4 (see page 111) 



126 


Download free eBooks at bookboon.com 



















Advanced Maths for Chemists 


Week 7: Chemistry and Vectors 


7 Week 7: Chemistry and Vectors 

We have been using vectors in Weeks 5 and 6 with matrices as the vertical column vectors (eigenvectors 
or wavefunctions). In our last two weeks we are going to use vectors more directly. In Fig. 7.1 the vector 
A is pointing in the positive directions along the x, y, z axes with components of a units, b units, c 
units, respectively. The vector is written as a column vector as below. The vector B is pointing in the 
negative x direction but the positive y and z directions which is shown in the signs of its column vector 
components. Some operations on vectors such as addition and subtraction may be treated purely as an 
exercise in trigonometry but we will use vector notation to link with the more complicated operations 
possible with vectors such as scalar and vector products. For clarity I will not explicitly show the unit 
vectors in the Cartesian coordinates. 
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7.1 Week 7 Tutorial Questions 

7.1.1 Questions 1: Time of Flight Mass Spectrometer 

A time-of-flight (TOF) mass spectrometer is sketched below. 


pulsed laser f 
sample 


A# 


detector 


V field-free flight tube 

Figure 7.2: schematic time of flight mass spectrometer. 


The TOF MS is all under a vacuum. A pulsed laser fires a pulse of light onto a sample plate causing a 
pulse of ions of mass m and charge q to be ejected from the sample (q = ze the number of charges z 
and e the elementary charge). The ions are accelerated by a voltage V applied between the sample plate 
and the flight tube. The flight tube is of length L (typically 1 metre) with no electric or magnetic fields 
in it. The pulse of ions leave the flight tube and strikes the detector. The electric field V gives the ions a 
kinetic energy E k with a velocity v of the ions in the flight tube. 

E k =^-=qV v = (2qV Imj 12 

The ions travel the distance L to the detector in a time t and rearranging the equations we have 

L 2 mll}\ m _ 2V 2 
v q\2 VI q L 2 

As the voltage V and the flight tube length L are fixed, the value of mlq is directly related to the square 
of the time of flight. Firstly, rewrite these equations in vector notation. Secondly, if L = 1 m, and V = 5 
kV, find the time of flight of a H + ion of mass 1.673xl0“ 27 kg and charge 1.602xl0“ 19 C down the flight 
tube (be careful with the units). The flight time in the short acceleration region may be safely ignored. 

Tump to Solution 1 (see page 131) 


7.1.2 Question 2: Dipole Moments of Molecules 

As vectors have magnitude and direction, the first thing to decide is the direction of the dipole. The 
electric dipole moment (or just dipole moment) (mu) is modelled as two equal but opposite charges 
in a molecule Q and -Q multiplied by a distance d separating the partial charges. 

A J= Qd 


128 


Download free eBooks at bookboon.com 












Advanced Maths for Chemists 


Week 7: Chemistry and Vectors 


It is a vector that points from the negative towards the positive end of the molecule. This is the internationally 
agreed sign convention. Some chemists use a non-standard arrow pointing towards the negative end to 
signify the movement of electrons and a cross at the positive end of the molecule. It should be made clear 
which system is being used, we will be using the “official” convention. The SI unit of dipole moment is 
the coulomb metre, which is much too large for use with molecules, so dipole moments are measured 
in debye (symbol D) where 1 D = 3.335641 xlO -30 C m. 


Cl 




Figure 7.3: chlorobenzene dipole moment. 


Chlorobenzene has a dipole moment of magnitude 1.5 D due to the polarized C-Cl bond. Calculate by vector 
addition the magnitudes and directions of the overall dipole moments for the isomers of dichlorobenzene 
orffio-dichlorobenzene (1,2-dichlorobenzene), meta-dichlorobenzene (1,3-dichlorobenzene) and para- 
dichlorobenzene (1,4-dichlorobenzene) assuming that the benzene ring is a perfect hexagon and the 
bond dipole [i remains constant at 1.5 D. 


lump to Solution 2 (see page 131) 


7.1.3 Question 3: The Absorption of Light by Polar Molecules 

It is possible to arrange molecules so that the bond dipoles all point in a particular direction. Fig.7.4 shows 
the alignment of carbonyl bonds in a sample. The magnitude of the carbonyl dipole moment is 2.5 D. 

y 

^ 31 ° 

x 


Figure 7.4: aligned carbonyl groups. 

The spectroscopic absorbance (A) of any sample is proportional to the magnitude of the component of the 
dipole which is parallel to the direction of polarization of the radiation used. Resolve the carbonyl dipole 
vector and hence calculate the ratio of the absorbances of the sample when the radiation is vertically 
polarized (v) to when it is horizontally polarized ( h ). 



lump to Solution 3 (see page 133) 
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7.1.4 Question 4: Hydrogen Halide Dipole Moments 

The magnitudes of the four hydrogen halide dipoles are given below. 



HF 

HC1 

HBr 

HI 

ft/ D 

1.91 

1.08 

0.80 

0.42 


Firstly using vector notation , express the magnitudes of the HF, HBr and the HI dipoles as scalar multiples 
of the HC1 dipole magnitude. Secondly Fig. 7.5 shows phosgene, in what direction does the phosgene 
molecular dipole point and what is its magnitude assuming the bond angles are 120° apart? 


Cl Cl 

\ / 


c 


II 


4= 1.5 D 



y B= 1.5 D 
30° 

C = 2.5 D 


Figure 7.5: phosgene. 


Tump to Solution 4 (see page 134) 
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7.2 Week 7 Tutorial Solutions 

7.2.1 Solution 1: Time of Flight Mass Spectrometer 


E k = 

2 

mv 

2 

--qV 

V = i 

2 qV 

V in 

L 

2 

m I 1 

m 

2 V 2 

V 

t = 

q 2 V 

q 

-u' 


1.673* 10' 27 kg _ 2*5* 10 3 V 2 
1.602* 10“ 19 C l 2 m 2 ’ 


1.673* 10~ 27 kgm 2 

( 1 . 602 * 10~ 19 cXl x 10 4 V) 


1.044* 10' 12 s 2 


t= 1.022* 10' 6 s 


Work = forcexdistance and soVC = J = Nm but also force = massxacceleration = kg m s 2 and so V 
C = J = kg m 2 s -2 . The proton has a flight time of about a microsecond. 


Return to Question 1 (see page 1281 


7.2.2 Solution 2: Dipole Moments of Molecules 



B = 


+ 1.5 cos 60 
— 1.5 sin 60 


\ 

+ 0.75 
-1.30 


Figure 7.6: 1,2-dichlorobenzene C-CI bond dipole vectors, units not shown. 


The two bond dipole moment vectors A and B each have their own components which is most clearly 
shown by drawing parallel Cartesian axes for the “dotted” B vector as shown in Fig. 7.6. Adding the two 
vectors gives us the molecules overall dipole vector C. 



'l.5] 


\ 

+ 0.75] 


\ 

+ 2.251 


C = 

0 

+ 

-1.30 

= 

-1.30 

|C| = V2.25 2 +( —1.30f+0 2 = 2.6 D 


0 J 


0 i 


0 i 



The angle, 0, that C makes with the v-axis is found from tan 0 = C y /C x as below, 


9 = tan 


-1.30D 
2.25 D 


= -30° 


131 


Download free eBooks at bookboon.com 

































Advanced Maths for Chemists 


Week 7: Chemistry and Vectors 


That is the molecules dipole moment bisects the bond dipoles as in Fig. 7.7. 



Figure 7.7: 1,2-dichlorobenzene dipole moment C = 2.6 D at 30° to both C-CI dipoles. 


The meto-dichlorobenzene vector diagram is shown in Fig. 7.8. 



Figure 7.8: 1,3-dichlorobenzene C-CI bond dipole vectors, units not shown. 


Adding the two vectors A and B gives us the molecules overall dipole vector C. 


C = 


/ 1.5 

0 

0 / 


-0.75 

-1.30 

0 


' + 0.75 1 
-1.30 


|C| = V0.75 2 + (-1.30) 2 + 0 2 = 1.5 D 


The angle, 9, that C makes with the x-axis is found from tan 9 = C /C x 


6 = tan 


- 1.30 D 
0.75 D 


= -60° 



Figure 7.9: 1,3-dichlorobenzene dipole moment C = 1.5 D at 60° to both C-CI dipoles. 


For para-dichlorobenzene the two bond dipole moment vectors are at 180° to one another and cancel 
out, adding the two vectors A and B gives us the molecules overall dipole vector C = 0 D. 
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1.5 


-1.5 

A = 

0 

B = 

0 


0 


0 


Figure 7.10: 1,4-dichlorobenzene C-CI bond dipole vectors, units not shown. 



' 1 . 5 ' 


'- 1 . 5 ] 


0 ) 

|c| = Vo 2 + 0 2 + 0 2 = OD 

c = 

0 

+ 

0 

= 

0 




0 ) 


1°) 



The dipole moment of 1,4-dichlorobenzene is 0 D, i.e. it is a non-polar molecule. 


Return to Question 2 (see page 128) 


7.2.3 Solution 3: The Absorption of Light by Polar Molecules 

In this question we are given the resultant molecular dipole moment vector and we have to resolve it 
into its horizontal ( h ) and vertical (v) component dipoles. 
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—2.5 sin 31 


-1.3' 

+ 2.5 cos31 

= 

+ 2.1 

0 i 


0 / 


Figure 7.11 : the components of the aligned carbonyl groups. 


The ratio of the spectroscopic absorbances for vertical polarisation (v) to horizontal polarisation ( h ) is 
independent of the signs of the components i.e. you cant have a negative absorbance. 


v _ |2.1 D| = 
h 11.3 Dl 


Return to Question 3 (see page 129) 


7.2.4 Solution 4: Hydrogen Halide Dipole Moments 

The ratios of the magnitudes are and thus the scalar multiples of the HC1 dipole magnitude are 


|A*hf| = 1.91 D , 77 
I^hci| 1.08 D 

I^HfI — 1*77 I^HCll 


1 ^HBrI _ 0.80D _ ~ 7z , 
I^hciI 1.08D 

I^HBrl = 0 * 74 I^HCll 


l^rnl _ 0.42P _ 
I^hci| 1.08 D 

| A*hi| — 0-39 |A^hci| 


The three bond dipole vectors of phosgene are shown below, 


A= 1.5 D 
60° > 


y B = 1.5 D 

+ 1.5 sin 60 


— 1.5 cos 30 


0 

II 

> 

O 

O 

CO 

\ 

X 

— 1.5 cos 60 

B = 

— 1.5 sin 30 

C = 

+2.5 

C = 2.5 D 

0 J 


0 i 


0 / 


Vector addition gives the overall dipole moment, D, of the molecule. 


A = 

' +1.3 ' 
-0.75 

B = 

' -1.3 ) 
-0.75 

C = 

0 ] 

+2.5 

thus D = 

0 

+ 1.0 


0 i 


0 i 


0 i 


0 / 


The horizontal components for each of the two C-Cl bond dipoles cancel one another out. The molecular 
dipole of phosgene is 1.0 D along the y-axis in the positive direction. 


Return to Question 4 (see page 130) 
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8 Week 8: Chemistry and Vectors 

8.1 Week 8 Tutorial Questions 

8 .1.1 Question 1: Methane Structure 



Figure 8.1 : methane in a cubic box of length 2 unit vectors. 
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Methane is tetrahedral, that is it fits inside a cube with the carbon atom at the intersection of any two 
diagonals, which is at the origin (0, 0, 0), Fig. 8.1. Two hydrogen atoms sit on the top face at the ends 
of a “horizontal” diagonal and the other two hydrogen atoms are on the bottom face at the ends of the 
other “horizontal” diagonal. Firstly, if the cube has sides of length 2 unit vectors calculate the bond 
vectors from the C-atom to each H-atom A, B , C and D. Secondly, choose any two bonds and use the 
vector product (scalar product) to calculate the angle between them. Thirdly, choose any two bonds and 
use the vector product (cross product) to calculate the area of the triangle formed by these two bonds. 

Tump to Solution 1 (see page 138) 


8 .1.2 Question 2: Magnetic Sector Mass Spectrometer 

A cation of charge q travelling with a velocity v enters a magnetic field B at right angles (9 = 90°) to the 
field and experiences a magnetic force F M which is the vector product, below. Draw a diagram to describe 
the direction of the force felt by a cation. What would be the direction of the force vector for an anion? 

q(v*B)= F m 

In a practical mass spectrometer the magnetic sector is not straight but the north and south poles are 
both curved round into an arc of a circle (typically 90°). The ion follows a circular path of radius r whilst 
inside this field. The centripetal force F c required to balance this magnetic force is F c = mv 2 /r. 

mv 2 i _ x 

-= q(v* B) 

r 

Draw a vector diagram of the various vectors and forces experienced by the positive ions during its time 
inside the circular magnetic field. 


Tump to Solution 2 (see page 140) 

8 .1.3 Question 3: Electron Spin and Orbital Angular Momentum 

From our atomic structure lectures we know an electron can only have has an intrinsic spin vector of s 
= V 2 . From quantum mechanics the magnitude of the spin vector is \s\ is, 

\s\ = V s(sr+1 ) h = V3/4 H 
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where h (h bar) is Planck’s constant divided by 2tt. Quantum mechanics also restricts the value of the 
component of s along the z-axis of the electron to m ^ to ± l Ah (spin up T and spin down >1). What about 
the other two Cartesian components of the spin vector (m x and m y ), can we say anything about them? 
Due to Heisenberg’s Uncertainty Principle if we know m z exactly then we can know nothing about m x 
and m y hence our spin vector points somewhere on a cone of uncertainty as in Fig. 8.2. Calculate the 
half-angle of the cone of uncertainty 9 of the electron spin angular momentum vectors. 



Figure 8.2: the two allowed orientations of the s = Vi electron spin angular momentum vector. 

The 2 p electrons have an orbital angular momentum l = lh. From quantum mechanics the magnitude 
of the orbital angular momentum vector is |/|. 

|/| = V7[7+Tj H = ^2H 

Quantum mechanics also restricts the value of the component along the z-axis to one of three possible 
values m z equal to +1 h, Oh , or - lh. Again because we know the value of the component along the z-axis 
exactly we know nothing about the values along the v-axis and y-axis. The orbital angular momentum 
vector points somewhere on one of three “cones of uncertainty” as in Fig. 8.3. Calculate the half-angles 
of the cones of uncertainty 9 of the electron orbital angular momentum vectors. 



Figure 8.3: the three allowed orientations of the / = 1 electron orbital angular momentum vector. 


Tump to Solution 3 (see page 140) 
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8.2 Week 8 Tutorial Solutions 


8 .2.1 Solution 1: Methane Structure 



Figure 8.4: methane C-H bond vectors. 


From the C-atom (at the origin) to the H-atoms the four bond vectors are given in Fig. 8.4. Secondly, if 
A and B are two vectors with components A x , A y A z and B x , B y B z with an angle 9 between them, then 
the scalar product (or dot product) may be written in either of two ways. 


A* B- A x B x +A y B y +A z B z = Ml|i?|cos0 The scaler product or dot product 


The scalar product as its name implies is a scaler quantity and is used most often in Chemistry to find 
the angle between two vectors and also for the potential energy of a magnetic dipole in a magnetic field 
( e.g . electron spin ESR and nuclear spin NMR). Choosing A and B as the two bond vectors (it doesn’t 


matter which pair you choose) we have for the left hand side of the scalar product. 



American online 

LIGS University 

is currently enrolling in the 
Interactive Online BBA, MBA, MSc, 
DBA and PhD programs: 


► 

► 


enroll by September 30th, 2014 and 
save up to 16% on the tuition! 
pay in 10 installments / 2 years 
Interactive Online education 

visit www.ligsuniversity.com to 
find out more! 


Note: LIGS University is not accredited by a 
nationally recognized accrediting agency li 
by the US Secretary of Education. 

More info here. 


138 



Download free eBooks at bookboon.com 























Advanced Maths for Chemists 


Week 8: Chemistry and Vectors 



-1 

+ i 

A* B = 

+ 1 

+ i 


r 1 / 

r 1 / 


(-lxl) + (lxl) + (-lxl) = -l 


For the right hand side of the scalar product we have 


^| = V(-l) 2 +l 2 + (-l) 2 = V3 |j?| = Vl 2 + 1 2 + 1 2 = V3 |4 14 cos6 = 3 cos6 


equating the two sides of the scaler product gives us the angle between the two C-H bonds in methane 
which is the tetrahedral angle. 


- 1 = 3 cos 6 0 = 109.5° 


Thirdly, we are asked for the area of the triangle formed by the C-atom and any two H-atoms, this is 
obtained by using the vector product or the cross product. An area can be represented by a vector at 
right angles to the surface of the parallelogram formed from the two vectors A and B as in Fig. 8.5. 




Figure 8.5: vector product or cross product and the right hand rule. 


,4X2* = W|S| sin d The vector product or cross product 


The direction of the vector AxB is found from the right hand rule. Arrange your thumb, index and 
middle fingers of your right hand at right angles to one another, align the thumb with the first vector 
A, then align the index finger with the second vector B , then the middle finger points in the direction 
of the third vector, the cross product |A||B|sin0, i.e. you count 1 , 2, 3 from your thumb for the three 
vectors A, B, and their vector product. 

AXB — V3xV3Xsin 109.5° = 2.828 units 2 


This the area of the parallelogram but we only want half of this for the triangle formed by vectors A and B. 


= 1.414 units 2 

2 
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8 .2.2 Solution 2: Magnetic Sector Mass Spectrometer 



thumb 


middl 



index 


Figure 8.6: direction of the magnetic force on a moving positive ion. 


q[vXB) — F m 


For a negative ion the force vector would point in the opposite direction. 



Figure 8.7: vectors and forces on a cation in a circular magnetic field. 


Return to Question 2 (see page 136) 


8.2.3 Solution 3: Electron Spin and Orbital Angular Momentum 
m = +14/? 


m= OH 

s 


m s = -!4 h 



= V 3 /i d 


Figure 8.8: the two allowed orientations of the s = Vi electron spin angular momentum vector. 





6 = cos 


1/2 a \ 
V 3 l 4 h) 


54.7° 


1 a -i 
m s = — — 6 = cos 


- 1/2 ft \ 

V3/4 n) 


125.3° 


The angles of the two half cones for the electron spin are measured from the positive z direction, so the 
m s = -Vi cone half-angle is 6= 180°-125.3° = 54.7° from the negative z- axis. The electron orbital angular 
momentum vector point in one of three orientations for the p atomic orbitals. 
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Week 8: Chemistry and Vectors 



Figure 8.9: the three allowed orientations of the / = 1 electron orbital angular momentum vector. 


m 


i 


= + 1 0 — cos -1 

l 


i n 

V2 H 


= 45° 


m l — 0 0 — cos 


0 n 
V 2 n 


= 90° 


rrij = — 1 6 = cos 


-i n 

V2 H 


= 135° 


The m { = -1 cone half-angle is 9 = 180°-135° = 45° from the negative z-axis. 

Return to Question 3 (see page 136) 

We have now finished our one year course of workbooks of Maths for Chemistry based students but 
keep using the three workbooks for revision and reference. Good luck with your degree. 



A cate-Lucent 
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In the near future, people may soon think it’s strange that 
devices ever had to be “plugged in.” To obtain that status, there 

needs to be “The Shift". 
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10 List of Formulae 


Advanced Maths for Chemists the following are some of the formulae encountered in this book. 


n 


ssr = Yjiyi ~ mx i~ c f 

i =1 

The sum of square of the residuals (SSR) for a straight line graph, 


Electrochemistry and thermodynamics, 

^ A H TAS 
nF nF 

Heat capacity of a solid using the Shomate equation, 

77 

C = A+Bt+Ct 2 + Dt 3 +—y 
p t 


Michaelis-Menten equation, 


Complex numbers, 


_ ^ max [S] 

^m + [S] 


i=(-l) L 


Euler s formula, 


exp(i@)= cos 6 + isin@ 


x' ] 


(1 

0 

0 ' 

/ 

X 

y' 

= 

0 

1 

0 

y 

z 'l 


1° 

0 

- 1 / 

\ z l 


Molecular symmetry operations, reflection in the xy-plane, 


( a — E 

\ p 


p 


\ 


a — E ) 


\ 

c i 



The Hiickel Matrix 
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a — E 

p 


p 

a—E 


= 0 


The Hiickel Secular Determinant 


A*B = A x B x +A y B y +A z B z = |A||B| cos@ The scaler product or dot product 
Ax B = | A|| B| sin@ The vector product or cross product 
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